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I.  INTRODUCTION 


This  is  the  final  progress  report  for  this  two  year  program  sponsored  by  the  Air  Force 
Office  of  Scientific  Research  on  "Domain  Processes  in  Ferroelectric  Ceramics".  As  an  issue 
which  has  both  scientific  and  technical  importance,  domain  process  in  ferroelectric  ceramics  has 
gain  more  and  more  attention  in  recent  years,  especially  in  the  area  of  material  development  for 
transducers  and  smart  structures.  People  are  eager  to  know  the  mechanism  of  domain  formation 
and  domain  dynamics  in  order  to  better  process  and  engineer  the  required  materials  to  achieve 
desired  physical  properties.  This  two  year  research  has  set  a  solid  foundation  to  build  upon  for 
future  studies  aimed  toward  more  thorough  understanding  of  domain  formation  and  domain 
dynamics. 

As  a  continuation  of  the  first  year's  effort  (Appendix  S>  we  have  made  substantial  new 
progress  in  the  second  year  as  described  in  this  report.  Due  to  the  complexity  of  the  problem  and 
limited  resources,  we  have  focused  on  a  few  key  issues  emphasizing  the  basic  theoretical 
development.  The  technical  progress  for  the  second  year  is  outlined  in  each  of  the  following 
sections,  and  the  details  can  be  found  in  the  four  appendices  A1-A4.  The  technical  report  for  the 
first  year  is  attached  as  Appendix  5  for  reference 

Section  n  describes  the  theory  development.  The  key  issue  is  the  nonlocal  coupling  of  the 
polarization  vector.  This  nonlocal  coupling  is  taken  into  account  through  a  set  of  polarization 
gradient  coefficients  in  the  energy  expansion,  which  eventually  enter  the  equilibrium  conditions 
in  the  Euler's  equation.  These  gradient  coefficients  are  proven  to  directly  couple  to  the  dispersion 
surface  of  the  soft  mode.  In  terms  of  the  microscopic  picture,  the  gradient  energy  is  the  mean 
field  representation  of  the  inter-site  coupling  of  electric  dipoles.  The  nonlocal  coupling  strength 
determines  the  domain  wall  width  and  the  coherent  length  in  domain  dynamics.  In  Section  in  we 
report  a  study  on  the  feasibility  of  using  electron  holography  to  ascertain  the  domain  profiles,  and 
try  to  extract  the  polarization  gradient  coefficients  from  the  electron  interference  pattern.  It  is  our 
belief  that  this  direction  is  worth  more  effort  in  die  future.  Section  FV  states  some  combined 
experimental  and  theoretical  effort  trying  to  identify  and  understand  more  complicated 
microdomain  structures  and  trying  to  find  out  the  intrinsic  mechanism  behind  the  modulated 
nanoscale  structures  in  doped  ferroelectrics.  Section  VI  highlights  the  current  development  and 
future  plans  which  are  generated  through  this  project,  and  are  still  carrying  on  the  momentum. 
Conqniter  simulation  on  the  domain  dynamics  is  beginning  to  take  off  and  we  believe  the 
dielectric  dispersion  in  the  domain  dynamics  will  be  understood  through  this  study. 

Many  results  were  produced  in  this  relatively  short  period  of  time  and  we  hope  the  seeds 
planted  here  will  flourish  in  the  near  future. 
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n.  DOMAIN  WALL  PROFILE,  POLARIZATION  GRADIENTS,  AND  THE  DISPERSION 

SURFACE  OF  SOFT  MODE 

Domain  walls  in  ferroelectrics  are  the  boundaries  of  different  variants  in  the  low 
temperature  ferroelectric  phase,  which  are  generated  during  the  paraelectric-ferroelectric  phase 
transition.  Due  to  the  long  range  nature  of  the  electric  and  elastic  interactions,  nonlocal  coupling 
of  the  electric  dipoles  is  very  strong,  and  this  nonlocal  coupling  determines  the  width  of  the 
domain  wails.  The  formation  of  a  spontaneous  polarization  can  be  well  characterized  by  a 
Landau-Devonshire  type  theory!  1].  Ferroelectric  systems  are  also  strongly  nonlinear  by 
including  the  nonlocal  coupling  in  the  free  energy  expansion  in  terms  of  polarization  gradients, 
one  can  derive  soliton-like  solutions  for  the  domain  walls[Al].  These  solutions  reflect  an 
intrinsic  mechanism  for  the  formation  of  domain  walls  without  interface  defects.  It  can  be  shown 
that  these  domain  states  have  higher  energy  than  single  domain  sate,  however,  they  can  be 
st^ilized  through  twinband  cross  pining  and/or  by  elastic  or  electric  boundary  conditions.  The 
formulation  described  in  Appendix  A1  has  included  the  elastic  coupling  as  well  as  the 
nonlinearity  and  nonlocal  coupling.  One  of  the  most  important  solutions  is  the  polarization  profile 
for  a  90°  domain  wall.  The  spontaneous  polarization  is  strongly  coupled  to  the  elastic  strain  in  a 
90°  domain  wall,  therefore,  is  directly  responsible  for  large  portion  of  the  field  induced 
piezoelectric  effect  in  ferroelectric  ceramics. 

Following  the  procedure  described  in  Al,  one  can  derive  the  temperature  dependence  of  the 
domain  profile  for  a  tetragonal  ferroelectric.  As  shown  in  Fig.  3  of  Appendix  Al,  with  the 
decrease  of  temperature,  the  magnitude  of  tlw  spontaneous  polarization  increases  but  the  domain 
wall  thickness  decreases.  The  changing  rate  fra*  both  the  spontaneous  polarizatirai  and  the  domain 
wall  thickness  are  most  noticeable  near  the  transition  temperature,  both  quantities  eventually 
become  insensitive  to  temperature  for  temperature  much  below  Tc.  Fig.  4  of  Appendix  Al 
depicts  the  influence  of  the  six  order  nonlinearity  parameter.  It  affects  very  strongly  on  the 
symmetrical  part  of  the  polarization  vector,  which  defines  the  local  structural  distortion  inside  the 
domain  wall ,  but  the  influence  on  the  antisymmetric  part  of  the  polarization  is  weak. 

The  most  influential  parameters  for  the  domain  wall  thickness  ate  the  polarization  gradient 
coefficients  which  reflect  the  strength  of  nonlocal  coupling.  This  effect  is  illustrated  in  Fig.  2  of 
Appendix  Al.  In  general,  stronger  nonlocal  coupling  will  have  broader  domain  wall.  This  is 
ctmceivable  since  the  domain  wall  is  actually  a  spatial  transition  region  between  two  coherently 
distorted  structures,  it  would  be  easier  to  transform  from  one  domain  state  to  the  other  when  the 
nonlocal  coupling  is  weak. 

Using  the  formulism  in  Appendix  Al  all  the  physical  quantities  associated  with  the 
ferroelectric  domain  walls  can  be  calculated.  However,  the  expansion  coefficients  in  the  free 
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energy  [Eqs.  (1-3)  of  Appendix  A  I]  must  be  determined  experimentally.  The  elastic,  dielectric 
and  electrostrictive  coefficients  can  be  directly  measured  and  are  readily  available  for  many 
materials.  The  difficulty  is  the  determination  of  the  polarization  gradient  coefficients,  which  many 
researches  only  estimate  through  some  hand  waving  arguments.  It  is  shown  in  Appendix  A2  that 
these  gradient  coefficients  can  be  directly  linked  to  the  microscopic  quantities  and,  in  principle, 
can  be  also  measured  directly.  A  simplified  rigid  ion  model  has  been  developed  for  a  perovskite 
system  as  described  in  Appendix  A2.  Since  the  polarization  which  is  the  density  of  electric 
dipoles  can  be  directly  related  to  the  relative  displacement  vector  of  positive  and  negative  charge 
centers,  the  gradient  coefficients  are  actually  directly  related  to  the  dispersion  surface  of  the 
corresponding  soft  mode  for  the  ferroelectric  transition.  Therefore,  through  inelastic  neutron 
scattering  experiments,  one  may  be  able  to  derive  these  polarization  gradient  coefficients.  This 
work  draws  a  clear  link  between  the  macroscopic  Landau  formulism  and  the  fundament 
microscopic  lattice  dynamics.  Although  the  problem  has  been  cleared  conceptually,  in  reality,  it  is 
far  from  completed  since  the  required  neutron  experiments  are  very  difficult  to  perform  due  to  the 
high  temperature  involved.  In  several  cases,  the  soft  mode  can  be  over  damped  which  makes  it 
impossible  to  measure.  A  recently  emerged  technique  Electron  Holography  offers  a  new  direction 
in  this  regard,  which  is  described  in  the  following  section. 

III.  EUECTRONHOLOGRAPHYANDTHEOBSERVATIONOFFERROELECrRIC 

DOMAIN  WALLS 

Electron  holography  utilizing  the  wave  characteristics  of  electrons.  Through  a  sharp 
emission  tip,  the  emitted  electron  beam  is  largely  coherent,  or  in  phase.  While  passing  through 
an  electric  field  region,  the  electron  wave  will  experience  a  phase  shift.  If  the  field  is 
inhomogeneous,  the  interference  pattern  formed  on  the  recording  film  will  be  twisted.  [2,3] 
Assuming  no  compensation  at  the  surface  of  the  electrodes,  a  twin  stracture  will  form  a  spatial 
varying  electric  field  with  large  field  gradient  in  the  domain  wall  region.  This  in  principle  will 
cause  a  twist  of  the  interference  pattern  produced  by  the  transmitted  electrons.  From  the  twisted 
interference  pattern,  domain  wall  thickness  c«i  be  evaluated  with  appropriate  data  interpretation. 
A  typical  twisted  interference  pattern  is  shown  in  Fig.  4  of  Appendix  A3,  which  is  produced  by 
coherent  electrons  transmitting  through  a  90'’  domain  wall  in  BaTi03.[3]  This  new  technique 
opens  up  a  different  avenue  to  obtain  the  gradient  coefficients  in  the  continuum  theory  described 
in  Appendix  Al,  since  one  can  do  back  fitting  using  the  experimental  data  and  the  theoretical 
derived  polarization  profiles.  However,  as  pointed  out  in  ref.  [4]  that  the  charge  compensation  at 
the  surface  of  the  ferroelectric  is  unavoidable,  therefore,  the  internal  field  supposed  to  be 
produced  by  the  spontaneous  polarization  is  actually  shielded,  so  that  the  interference  pattern 
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twisting  will  not  be  seen  in  many  cases.  These  twisted  interference  patterns  can  be  seen  only  in 
some  situations  when  the  compensation  time  constant  is  very  long,  or  momentarily  with  the 
disturbance  of  temperature.  [4]  More  thorough  theoretical  understanding  is  still  needed  to  perfect 
the  interpretation  of  the  electron  interference  patterns. 

Some  charged  defects  inside  the  sample  can  also  produce  disturbance  to  the  interference 
pattern  through  the  electric  field  surrounding  these  defects.  Fig.  5  in  Appendix  A3  shows  this 
situation.  An  interesting  point  is  the  interaction  of  these  charged  defects  with  domain  walls,  it 
s^pears  that  they  are  attracted  to  the  domain  wall  region.  This  observation  provides  an  important 
evidence  for  the  defect  pinning  to  the  domain  walls,  which  gives  us  a  possible  explanation  of 
why  a  few  percent  of  aliovalent  doping  could  change  the  physical  properties  of  piezoelectric 
ceramic  PZT  so  drastically[A3].  Many  potential  applications  of  this  new  technique  await  for 
further  exploration. 

IV.  MICROSTRUCTURAL  MODULATIONS  IN  PLT 

In  domain  dynamics,  the  inertia  associated  with  domain  wall  motion  is  quite  large  since  the 
whole  domain  has  to  be  moved  in  the  perpendicular  direction  to  the  domain  wall  motion.  It  is 
conceivable  that  the  domain  wall  motion  will  become  easier  when  the  domain  size  becomes 
smaller.  Small  domains  may  be  achieved  by  aliovalent  doping  since  the  charged  dopants  have 
strong  interaction  with  domain  walls,  hence  will  affect  the  domain  formation.  Smaller  domain 
size  makes  it  easier  to  switch  and  produces  relatively  larger  dielectric  response.  Appendix  A4 
reports  a  study  on  Lanthanum  doped  lead  titanate  which  produces  a  spectrum  from  pure 
ferroelectric  to  relaxor  ferroelectric.  As  shown  in  Fig.  3  of  Appendix  A4,  the  domain  structure 
changes  from  regular  large  size  at  1%  lanthanum  doping  to  an  ill-defined  mirostructural 
modulation  at  25%  doping.  The  strongly  first  order  ferroelectric  phase  transition  in  PbTi03  is 
also  gradually  changed  to  a  nearly  second  order  transition.  The  emerge  of  microdomain-like 
structures  within  a  well  defuied  domain  during  the  increase  of  lanthanum  doping  provides  a  key 
for  the  undostanding  of  dopant  driven  transitimi  from  a  conventional  ferroelectric  to  a  relaxor 
ferroelectric.  The  domain  signature  gives  us  a  direct  way  to  distinguish  between  different  types 
of  ferroelectrics,  the  dielectric  and  piezoelectric  performance  of  ferroelectric  materials  are  also 
directly  coupled  to  the  domain  structures. 
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VI.  PROGRAMS  IN  PROGRESS 

Although  the  cuirent  two  years  contract  has  come  to  a  conclusion,  the  effort  initiated  by  this 
Air  Force  grant  has  planted  a  seed  for  the  future  study  of  domain  processes  in  ferroelectrics.  The 
understanding  on  this  subject  is  crucial  in  the  development  of  functional  ceramics.  One  of  the 
natural  extensions  of  this  work  is  to  study  domain  dynamics.  According  to  our  analyses,  the 
problem  may  be  simplified  to  a  two-dimensional  electric  dipole  problem  with  inter-site  coupling. 
Therefore,  molecular  dynamics  methods  can  be  used.  We  are  currently  developing  a  computer 
program  based  on  the  formulation  described  in  Appendix  A I  and  adding  an  alternating  electric 
field.  At  the  beginning  stage,  only  one  twin  will  be  calculated,  and  later  we  plan  to  extend  the 
calculation  to  multidomain  configurations.  The  primary  focus  is  on  the  90°  twin  for  which  the 
elastic  strain  is  directly  coupled  to  the  domain  wall  motion.  It  is  anticipated  that  this  calculation 
can  shed  light  on  the  extrinsic  contributions  to  the  dielectric,  piezoelectric  and  electrostrictive 
effects  in  ferroelectric  ceramics. 

We  also  plan  to  carry  out  more  experimental  works  using  electron  holography  technique, 
which  not  only  helps  to  develop  a  powerful  new  technique  for  the  study  of  domain  ferroelectrics, 
but  also  can  push  forward  the  theoretical  study  on  domain  walls  as  well  as  understanding  the 
interaction  of  domain  walls  and  aliovalent  dopants. 

Many  ambiguous  points  were  cleared  through  this  study,  such  as  the  morphotropic  phase 
boundary  in  P2T,  the  development  of  microdomatn  structures  through  aliovalent  doping,  and  the 
relationship  between  the  macroscopic  Landau  theory  and  the  microscopic  lattice  dynamics.  A 
new  technique  for  direct  observation  of  the  domain  wall  profile  is  analyzed  and  some  theoretical 
ground  works  are  set  for  the  proper  interpretation  of  electron  interference  patterns  generated  by 
coherent  electron  going  through  ferroelectric  twins. 

Many  problems  still  exist  in  the  domain  studies,  and  it  is  our  belief  that  the  understanding 
of  these  domain  processes  not  only  will  enhance  our  knowledge  about  this  fascinating  mateiial 
structure,  but  also  can  provide  firm  theoretical  background  for  the  future  development  of  better 
functional  ceramics. 
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Abstract  The  domain  structures  in  ferroelectrics  can  be  described  by  a 
Landau-Ginzburg  type  theory  with  the  twin  and  twin  band  (domain)  structures 
being  nonlinear  and  nonlocal  excitations  of  the  ferroelectric  phase.  The 
polarization  gradients  in  the  theory  reflect  the  degree  of  nonlocal  coupling 
along  different  crystallogr^hic  orientations.  These  gradient  parameters  can  be 
obtained  either  from  the  dispersion  surface  of  the  soft  mode  or  through  fitting 
the  polarization  profile  measured  by  the  holographic  electron  microscopy. 


INTRODUCTION 

The  understanding  of  domain  structures  is  essential  for  the  design  and  applications  of 
ferroelectrics.  It  has  been  long  recognized  that  the  piezoelectric  and  dielectric  properties 
of  ferroelectric  ceramics  are  mainly  determined  by  the  behavior  of  domain  structures. 
The  formation  of  domains  in  ferroelectrics  is  due  to  the  existence  of  milti-variants  in  the 
ferroelectric  phase.  Atomic  coherency  is  usually  maintained  across  the  domain 
boundaries,  which  make  it  possible  to  switch  domain  orientations  from  one  to  the  other 
using  external  (either  mechanical  or  electrical)  fields.  This  switching  gives  rise  to  the  so 
called  extrinsic  contributions  to  the  materials  properties.  The  formation  of  domain  walls 
in  ferroelectrics  may  be  treated  in  terms  of  solitory  wave  excitations  in  a  nonlinear  and 
nonlocal  system.  Single  kink-like  and  periodic  solitory  wave  solutions  for  the  twin  and 
periodic  domain  structures  can  be  derived  using  the  continuum  theory.  Since  all  the 
expansion  coefficients  in  the  Landau-Devonshire  model  can  be  expressed  in  terms  of 
measurable  macroscopic  quantities,  the  continuum  theory  can  give  quantitative 
description  of  die  domain  wall  properties,  including  the  profile  of  polarization  across  the 
domain  wall,  domain  wall  width,  energy  stored  in  the  muti-domain  structure,  and  the 
stress  build  up  at  the  domain  wall  region,  once  the  polarization  gradient  coefficients  are 
obtained. 


THE  MODEL 


The  Landau-Devonshire  type  phenomenological  theory  for  ferroelectrics  has  been 
developed  for  the  ferroelectric  phase  transition. For  a  cubic  system,  such  as 
perovskite  ferroelectrics,  the  elastic  Gibbs  free  energy  can  be  expressed  in  the  following 
form: 

GsGp  +  Gel  +  Gc  (1) 

Gp  =  A  (Pi2  +  P22  +  P32  )+  B  (Pi^  +  P2'‘  +  +  C  (P16  +  P26  +  P36  ) 

+  D  (Pi2  P22  +  P22  P32  +  P,2  P32)  +  E  (  Pi^  P22  +  P,2  P24  +  P24  P32  +P22  P34 
+  P32p,4  +  p,2p34)  +  HPi2p22p32  (2) 

Gel  =  -  ^  ( Xi  i2  +  X222  +X332 )  -  S12  (  XnX22  +  X22X33  +  X,  1X33) 

-"-f  (X122  +  X132+X232)  (3) 

Gc  =  Qi  1  ( Xi  1  Pi2  +X22  P22  +X33  P32 )  +  Qi2  [  Xn  (P22  +  P32 )  +  X22  (Pi2  + 
P32  )  +  X33  (Pl2  +  P22  )1  +  Q44  (  X12  Pi  P2  +  Xi3  Pi  P3  +  X23  P2  P3  ) 

(4) 

where  A,  B,  C,  D,  E,  H  are  the  linear  and  nonlinear  dielectric  constants,  sij  are  the 
elastic  compliance  coefficients,  Qij  are  the  electrostriction  constants.  Pi  and  Xij  are  the 
components  of  polarization  and  stress,  rer.pectively.  All  the  coefficients  are  assumed  to 
be  independent  of  temperature  except  A  which  is  a  Uneariy  function  of  T, 

A  =  a(T.To)  (5) 

In  a  homogeneous  system,  a  paraelectnc-ferroelectric  phase  transition  occurs  at  Tc. 
Under  stress  free  condition,  the  phase  tra’^sition  temperature  Tc  and  the  spontaneous 
polarization  Pc  at  the  transition  can  be  deriv::^  by  mi  >imizing  Eq.  ( 1 ), 

Tc  *To  +  -^  (6) 

4  C  a 

Pc^-^.  (7) 

One  of  die  low  temperature  ferroelectric  phases  is  the  tetragonal  phase.  There  are  six 
energetically  degenerate  vahvints  in  the  tetragonal  phase:  (±  Ps,  0, 0 ),  (0,  ±  Ps,  0 )  and 
(0, 0,  ±  Ps),  where  P*  is  rhe  spontaneous  polarization  given  by 

<«> 

These  energetically  degenerate  variants  can  coexist  in  the  ferroelectric  phase  to  form  the 
twin  strucft’.res.  Electron  microscopy  reveals  that  the  ionic  coherency  is  maintained 
across  domain  walls,  but  domain  walls  are  not  atomically  sharp.  Domain  wall  width  is 
determined  by  the  nmilocal  coupling  strength  of  the  ferroelectric  system. 


Since  the  nonlinearity  has  been  included  in  the  model  [see  Eqs.  (l)-(4)],  if  we  add 
the  contributions  of  nonlocal  coupling,  then  from  soUton  theory,  we  may  expect  solitory 
wave  excitations  in  the  ferroelectric  phase.  These  excitations  are  in  deed  found  and  they 
represent  the  ferroelectric  domain  walls. 

For  a  perovskite  system,  the  symmetry  of  the  high  temperature  phase  is  cubic, 
therefore,  the  Gibbs  energy  representing  the  nonlocal  coupling  can  be  written  as 
follows:^ 

Gg  =  5gll  (Pl.l^  +  P2,2^'^P3,3^)  +  gl2(Pl.l  P2.2'^Pl.l  P3.3 P2.2  P3,3  ) 

+  ^  g44  [  (Pl.2  +  P2.i)2  +  (Pl.3  +  P3.l)^  +  (P2.3  +  P3.2)2]  (9) 

here  the  indices  after  the  comma  represent  derivatives  with  respect  to  space  variable 
along  that  axis.  Upon  the  minimization  of  the  total  free  energy  of  the  system  Eqs.  ( 1) 
and  (9),  one  can  obtain  the  solutions  for  the  domain  walls.^ 


90°  DOMAIN  WALLS 

There  are  two  types  of  domain  walls  in  the  tetragonal  ferroelectrics.  One  is  the  180° 
domain  wall  which  divides  a  twin  domain  with  identical  strain  but  opposite  polarization, 
and  the  other  is  the  90°  domain  wall  which  divides  two  domruns  whose  polarization  and 
spontaneous  strain  ate  nearly  90°  from  each  other.  Solutions  for  the  180°  domain  walls 
can  be  easily  obtained  since  the  problem  is  one-dimensional.^  Here  we  only  solve  the 
problem  of  90°  domain  walls  for  which  the  problem  can  be  rendered  to  quasi-one- 
dimensional. 

From  transmission  electron  microscopy  studies,  domain  walls  tend  to  broaden  or 
bent  near  the  surface,  however,  inside  the  sample  they  have  well  defined  wall-like 


FIGURE  1.  A  tetragonal  twin  structure  and  the  coordinate  system  used  in  this  paper. 


structure  with  translational  symmetry  parallel  to  the  wall  plane.  Therefore,  while  dealing 
with  a  <llO>-type  domain  walls  in  tetragonal  ferroelectrics  we  can  rotate  the  xi,  X2 
coordinates  about  the  X3  coordinate  by  45°  so  that  the  properties  of  the  domain  walls 
only  depend  on  one  space  variable  (s-coordinate  as  indicated  in  Fig.  1)  only. 


In  the  new  coordinate  system  the  equilibrium  conditions  are  governed  by  the 
following  equations:^ 


8  8G  30 

=0,  (i.j  =  s,r,3) 

3xj  3Pi,j  3Pi 

(10) 

Xjjj  =  0,  (i,j  =  s,  r,  3) 

(11) 

and  we  also  need  the  elastic  compatibility  rel^ons 

Eikl  £jmn  xin,iun  “  ®  (i,j)k,l,m,n  =  s,  r,  3) 

(12) 

to  insure  the  elastic  continuity  since  in  our  model  the  domain  walls  are  intrinsic 
excitations,  no  defects  are  created  in  the  domain  wall  region.  xin  is  the  component  of 
elastic  strain  tensor  and  eiu  is  the  Levi*Civita  density. 

Eqs  (1 1)  and  (12)  has  three  nontrivial  solutions: 

Xr3  =  0  (13) 

Xn* - ^ 2S12Q12  *Sii(Qii+Qi2)1Po‘1  2S12Q12  -Su(Qu+Ql2-Q44)lPs 

2  (sns„-st2) 

-[2S12Q12 -sn(Qn+Qi2+Q44)lPr }  (14) 


X33= - ^ - - — {[S12(Qii+Qi2)  -2S5sQi2]Po-[  Si2(Qi  I+Q12+Q44)  -2SssQl2lPr 

2(si  1853-8]  2) 

-[Sl2(Qn+Ql2-Q44)  -2SssQi2]  P5}  (15) 

where 

Ss5*^S11+S12+^)- 

It  can  be  easily  verified  that  the  two  stress  components  Xrr  and  X33  are  nonzero  only  in 
the  vicinity  of  the  domain  wall.  These  nonzero  stress  components  near  the  domain  wall 
region  is  the  cause  of  the  faster  etching  rate  which  makes  the  domain  walls  visible 
through  chemical  etching  technique. 

In  order  to  see  the  general  trend  of  the  variation  of  polarization  proHle  without 
specifying  the  coefficients  to  a  particular  system,  we  normalize  the  polariztUion  and  the 
space  variable  s  into  dimensionless  forms  by  the  following  substitutions: 


say§. 


=  P.=f^  f.  =  Pcf, 


where 


(16a,b) 

(I7a,b) 


G$s*^gU+gl2+2g44).  Grs=^gn-gl2). 
and  define  the  dimensionless  temperature  as 

^Tc-To 


(18) 


then  the  equilibrium  condition  Eq.  (10)  can  be  written  in  the  following  form  for  a  90° 
twin  structure. 


a  f,,  Xs  fs+  bsf|+  c  fsf?+  dff  +  (8  -  ftf?+  (4  -  ^)  f^f? 
i  f,.  Xr  fr+  brf?+  C  frff  +  df?+  (8  -  ^)  f3f|+  (4  -  ^)  f.ff 


where  the  coefficients  are  given  by 


-V5 

(l+i/  I- it) 

Cs=t - * - J - 

3(SiiSss*St2)B 


{(Q1I+Q12  -  Q44)I2S|2Qi2-Sii(Qii+Qi2)1 
+2Qi2[S12(Qi1+Qi2)-2S5,Qi2]} 


x^- 


jl+^l  l-la) 

3(siiS„-sf2)B 


{  (Ql  1+Qi2+Q44)[2si2Qi2-Si  i(Qi  1+Qi2)1 
+2Qi2[S12(Qi  I+Qi2)'2s5,Qi2]  } 


(19) 

(20) 


(21) 


(22) 


(23) 


b.>-2-D4 


-{  (Ql  l+Ql2-Q44)[2Si2Ql2-Sl  l(Ql  1+Qi2-Q44)1 


B  2B(snSss-sf2) 

+2Qi2[Si2(Ql  1+Qi2-Q44)-2s„Qi2]  } 


(24) 


b,.-2-P.t-— .-1 


-{(Qu+Qi2+Q44)[2Si2Qi2-Si  l(Ql  I+Q12+Q44)] 


B  2B(sns„-sf2) 

+2Qi2[Si2(Qi  l+Qi2+Q44)-2SssQi2]  } 


(25) 


± 


-{(Qu+  Q12-  Q44)[2Si2Qi2-  S11(Qi1+Qi2+Q44)1 


B  2B(sus„-sy 

+2Qi2(Si2(Qll+Ql2+Q44)-2s$$Qi2l } 


(26) 


(27) 


Using  the  normalized  equations,  we  can  study  the  influence  of  different  parameters  to 
the  polarization  profile  and  obtain  a  conceptual  understanding  on  the  nature  of  the 
polarization  variation  in  the  domain  wall  region.  Fig.  2  shows  the  variation  of  the 
polarization  components  with  the  parameter  a  across  a  90’’  domain  wall.  One  can  see 
that  the  domain  wall  becomes  wider  as  a  increases.  In  real  dimensions,  because  the 
scaling  factor  of  the  space  variables,  %  is  directly  related  to  the  product  Gss*Grs  [see 
eq.(I7)],  domain  wall  becomes  wider  as  the  gradient  coefficients  become  larger. 


FIGURE  2.  Variation  of  polarization  components  fs  and  fr  induced  by  the  change  of 
parameter  a  across  a  90°  (k»nain  waU.  The  gradient  parameter  a  determines  the  width 
of  the  domain  wall. 


FIGURE  3.  Variation  of  polarization  c  mponents  fs  and  fr  with  temperature  t  across 
a  90°  domain  wall.  The  asymptotic  values  of  die  polarization  components  increase 
and  the  domain  wail  width  decrease  while  lowering  temperature. 
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FIGURE  4i  Variation  of  polarization  components  fs  and  fr  induced  by  the  change  of 
parameter  d  across  a  90**  domain  wall.  The  nonlinear  parameter  d  influences  the 
magnitude  of  the  polarization  variation  in  the  domain  wall  region. 
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FIGURE  5.  (a)  Polarizatitm  components  fi  and  f2  across  a  90"  domain  wall,  and  (b) 
Ohistration  of  die  variations  of  the  polarization  vector  and  the  unite  cell  distention 
across  a  90"  domain  wail. 


Fig.  3  shows  the  variation  of  the  polarization  with  temperature  t.  The  asymptotic 
values  of  the  magnitude  of  the  polarization  components  increase  and  the  domain  wall 
thickness  decreases  as  the  temperature  is  lowered.  The  temperature  dependence  is  strong 
near  the  transition  and  gradually  becomes  insensitive  when  the  temperature  is  far  below 
Tc.  Fig.  4  shows  the  variation  of  the  polarization  components  fs  and  fr  induced  by  the 
change  of  parameter  d.  We  can  see  that  the  magnitude  of  fs  is  very  sensitive  to  d  while 
the  domain  wall  thickness  is  relative  insensitive  to  d.  We  can  also  calculate  the 
polarization  components  fi  and  (2  in  the  original  coordinates.  One  example  is  given  in 
Fig.  S  (a)  for  a  set  of  chosen  parameters.  The  corresponding  unit  cell  distortion  and  the 
polarization  variation  across  the  domain  wall  are  illustrated  in  Fig.  5(b).  The  polarization 
vector  rotates  gradually  from  one  orientation  into  the  other  accompanied  also  by  a 
change  of  the  magnitude. 

All  properties  of  the  domain  walls  can  be  quantitatively  calculated  using  this 
model  (Mice  the  expansion  coefficients  are  known.  As  we  have  mentioned  above  that  the 
polarization  gradient  coefficients  are  most  crucial  for  the  study  of  domain  walls,  which 
may  be  derived  from  the  measurements  on  the  dispersion  surface  of  the  soft  inode.^>^  In 
general,  inelastic  neutron  scattering  to  probe  the  soft  mode  may  be  difficult  due  to  the 
relatively  high  transition  temperature  in  many  systems  of  interest  and  in  some  cases,  the 
soft  mode  is  over  damped.  An  alternative  way  to  obtain  these  coefficients  would  be  to 
probe  the  polarization  profile  across  the  domain  wall  and  then  fitting  the  unknown 
coefficients  using  the  differential  equaions  (19)  and  (20).  The  recently  emerged  new 
technique,  electron  holography,  may  offer  an  (^ti(Mi  to  this  end.^^'^  > 
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Abstract 

The  gradient  coefficients  in  the  Landau-Ginzburg  theory  are  crucial  for  quantitative 
description  of  domain  walls  in  ferroelectrics.  The  magnitude  of  these  gradient  coefficients  are  a 
measure  of  nonlocal  coupling  strength  of  the  polarization.  In  this  paper,  we  intend  to  explain  the 
physical  meaning  of  these  gradient  coefficients  in  terms  of  lattice  dynamics  and  give  some 
relaciofiships  between  these  gradient  coefficients  and  the  dispersion  surface  of  the  soft  mode.  The 
implications  for  the  study  of  over  damped  soft  modes  are  also  discussed. 

PACS  Numbers:  64.60.-i,  63.20.Dj,  64.90.+b 
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I.  Introduction 

Many  fem)electnc  materials  have  perovskite  structure  with  a  cubic  symmetry  in  the 
paraelectric  phase.  The  symmetry  of  the  low  temperature  ferroelectric  phase  can  be  tetragonal, 
rhombohedral  or  orthorhombic  Ferroelectric  phases  usually  have  more  than  one  variants  and 
these  variants  may  coherently  coexist  within  the  synunetry  frame  of  the  parent  phase,  forming  the 
so  called  twin  structures.  It  is  shown  that  these  twin  structures  can  be  well  described  by  Landau- 
Ginzburg  (LG)  type  models.  L2  All  the  expansion  coefficients  in  the  Landau  theory  correspond  to 
certain  macroscopic  physical  quantities  and  can  be  obtained  experimentally.  However,  the 
physical  meaning  of  the  gradient  coefficients,  which  regulate  the  domain  wall  formation  and 
control  the  domain  wail  width  in  the  twin  structures,  still  needs  to  be  specified. 

A  paraelectric-ferroelectric  phase  transition  is  characterized  by  a  softening  of  a  transverse 
optic  mode  at  the  Brillouin  zone  center  due  to  the  cancelation  of  the  long  range  Coulomb  forces 
and  the  short  range  repulsive  forces.^  The  soft  mode  is  stabilized  above  the  phase  transition 
temperature  by  the  anhaimonic  interactions  whose  strength  weakens  as  the  temperature  decreases. 
Using  mean  tield  theory,  one  can  still  formally  retain  the  terminology  of  normal  nxides  if  the 
"soft  mode"  frequency  is  assigned  to  be  temperature  dependent. 

Qose  to  the  phase  transition  temperature,  the  dominant  contribution  in  the  lattice 
Hamitimian  is  from  the  soft  mode.  Therefore,  one  may  simply  study  the  soft  mode  behavior  to 
characterize  the  phase  transition  near  Tc.  For  an  inhomogeneous  system,  the  spatial  variation  of 
the  order  parameter  must  be  considered.  This  is  done  by  adding  an  energy  term  induced  by  the 
order  paramrter  gradients.  Because  the  inclusion  of  both  nonlinear  and  nonlocal  terms  in  the 
energy  expansion,  oac  may  expect  to  obtain  large  amplitude  soliton-like  solutions  which  can 
describe  the  domain  walls  (the  transition  region  between  coherent  twin  structutes).^'^’^  The 
physical  nr««»ning  of  these  gradient  coefficients  can  be  seen  from  the  study  of  small  amplitude 
oscillations  for  which  the  nonlocal  coupling  can  be  treated  as  perturbations.^*'^  We  will  show  in 
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this  paper  how  these  polarization  gradients  can  be  derived  from  a  simplified  lattice  dynamical 
model. 

Since  polarization  is  the  density  of  dipoles  per  unit  volume,  it  is  proportional  to  the 

magnitude  of  the  associated  optical  mode.  As  will  be  shown  in  the  third  section  of  this  p£^r  that 

the  lattice  potential  for  a  given  optical  mode  can  be  written  in  terms  of  the  polarization  vector. 

Therefore,  the  Landau-Ginzburg  potential  can  be  directly  used  in  the  lattice  dynamical  calculations 

in  the  small  representation  of  the  soft  mode. 

n.  Gradient  Coefficients  and  Dispersion  Surface  of  the  Soft  Mode 

From  soft  mode  theory,  the  potential  energy  for  a  cubic  system  may  be  expanded  in  terms 

of  the  eigenvector  of  the  soft  mode^ 

G  *  ^  (uf +u^ +u^  )  (1) 

where  k  is  related  to  the  temperature  dependent  soft  mode  frequency,  k  oe  (T-Tc),  and  Ui  (i  = 

1 ,2,3)  are  the  con^nents  of  the  eigenvector  of  the  soft  mode.  For  the  ferroeiectric  phase 
transition  u  is  a  relative  displacement  field. 

If  tt  is  inhomogeneous,  we  must  include  die  gradient  energy  in  the  energy  expansion.  For 
cubic  symmetry  the  gradient  energy  may  be  written  as  follows, 

Gg  =  ^  (u^.l  +u5.2  +«5.3  )  +Sl2  (Ul.l  U2.2  +  Ui,i  U3.3  +  U2.2  U3.3  ) 

(  (  Ui,2  +U2.i)2  +  (  Ui ,3  +U3.i)2  +(  U2.3  +U3.2)^  )  (2) 

Assuming  the  effective  mass  for  the  mode  is  M  then  the  equations  of  motion  become 
M  tti  +  K  UI  -  811  Ui.il  -  812  (U2,21  +  U3,31  )  -  844  (Ul,22  +U2,12  +  Ul,33  +  U3,13  )  =  0,  (3a) 

M  U2  +  K  U2  -  811  U2,22  -  812  (Ul,12  +  U3.32  )  -  844  (Ul.21  +U2,H  +  U2,33  +  U3.23  )  *  0.  (3b) 

M  tt3  -f  K  U3  -  81 1  U3,33  -  812  (Ul,i3  +  U2,23  )  -  §44  (Ul,31  +U3,n  +  U2,32  +  U3,22)  =  0.  (3c) 

Eq.(3a<)  have  plane  wave  solution  of  the  fwm 

tt  s  U  exp  [  j  (OK  -k*x)l  (4) 

Substituting  eq.  (4)  into  eqs.  (3a-c)  gives  the  eigenvalue  problem 
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D(k)  U 

where  D(k)  is  the  dynamical  matrix 


(5) 


D(kH 


K+5iik|+644(k2+k3)  612  ki  k2 

812  ki  k2  K+6iiki+844(k?+k3) 

8i2  ki  k3  612  k2  k3 


612  ki  k3 
5i2  k2  k3 


(6) 


ic+6iikf+544(k?+k|)  / 


If  the  depolarization  field  is  included,  the  equations  of  motion  (3a-c)  will  contain  one  more 
term  representing  this  contribution,  which  will  split  the  longitudinal  and  transverse  optical 
inodes.^*^  The  depolarization  field  is  given  by 


E(k)=r 


(P*k)  ^ 

Co  k^ 


(7) 


The  additional  contribution  is  a  linear  function  of  the  polarization  vector  P  which  is  proportional 
to  the  relative  displacement  field  u.  Adding  eq.  (7)  to  the  r.h.s.  of  eq.  (3a-c)  leads  to  the 
dynamical  matrix  for  a  given  k  hi  what  follows,  we  will  treat  three  k>values  in  the  three  principle 
directions  of  the  k>space. 

A.  ks[k.O.O]. 

Define  PisZ  e  Ui,  where  Z  is  a  constant  which  has  a  unit  of  inverse  volume  and  e  is  the 
electron  charge  unit  The  meanigng  of  Z  will  be  clear  fiom  later  derivations.  For  this  k  value,  the 
dynamical  matrix  can  be  sin^riified  to  the  following  form 


K+A  +  8iik^ 

0 

0 

D(k)* 

0 

K+844  k^ 

0 

0 

0 

K+844  k^ 

where  A  is  a  constant  defined  by 

A  *(l/eo)  (Ze)2 

(9) 
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From  Eq.  (8)  one  can  easily  derive  the  dispersion  relations  for  the  longitudinal  (qjl)  and 
transverse  ((Of)  modes  respectively: 


ci)f  s^(K  + A  +  811  k2)  (10a) 

Q)t  =  jJi  (K  +  844  k2  )  (10b) 


B. 


k 

VI 


[1.  l.Ol 


The  electrostatic  force  from  the  depolarization  field  is  now  given  by 


-  -^^2)  [11^  0]  a  ^  “2)  [1. 1. 0] 

2  eo  ^ 

Therefore  the  dynamical  matrix  is 


(11) 


D(k)s| 


K+A/2+  (811+844)  k2/2 
A/2  +(8i2/2) 

0 


A/2  +  812  k^ 
K+A/2+  (811+844)  k2/2 
0 


and  the  dispersion  relations  are 


0 

0 

K+  844  k^ 


(12) 


<aj  = 

=  jL(K  +  A  +  ^(8li+8i2  +  844)k2] 

(13a) 

4, 

*  ^  t  K  +  ^  (811  -812  +  844 )  k^  ] 

(13b) 

4 

=  iJi[K+  844k2] 

(I3c) 

Here 

the  two  transverse  modes  are  not  degenerate. 

C.  k 

[1. 1.  n. 

Fot  this  case  the  Coulomb  force  from  the  (^polarization  field  is 

Z  e  (PLt  P2-*-  P3}  [I  ,A  (Ui+  U2+U3)  ^ 

3  eo  ^ 

and  the  dynamical  matrix  becomes 
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Dik)4 


k+A/3  "KSi  1/3+2  644/3) 
A/3  +  (612/3)  k2 
A/3 +  (612/3)  k2 


A/3 +  (612/3)  k2 
k+A/3  +(611/3+2  644/3)  k’ 
A/3  +  (612/3)  k^ 


A/3  +  (612/3)  k- 
A/3  +  (612/3)  k‘ 
k+A/3 +(611/3+2  644/3)  k- 


(15) 


The  dispersion  relations  are  therefore  given  by 

ai  =X[K  +  A  +  l(6u+2  6,2+2644)k2]  (16a) 

<•4  3(^11 -812 +  2644)  k^  ]  (16b) 

Note  that  the  dispersion  relations  derived  £dx>ve  are  for  the  cubic  phase  near  k=0  but  not  for 
the  low  temperature  ferroelectric  phase. 


m.  The  Expansion  Coefficients  and  Lattice  Dynamics 

Taking  the  limit  k  0  in  the  dispersion  relations  derived  above,  one  finds  that  the 
coefficient  k  /  M  becomes  the  soft  transverse  mode  frequency  square, 

=  k/M  (17) 

k^O  ^ 

The  longitudinal  mode  will  not  become  soft  at  T  s  Tc  because  of  the  depolarization  field 

contribution  A, 


^0  <'8) 

The  sinqslest  model  to  calculate  these  coefficients  in  terms  of  microscopic  quantities  is  to 
study  the  k  aO  mode  for  a  biatomic  system  using  rigid  ion  model,  in  which  the  soft  mode 
represents  the  relative  displacement  field,  the  mass  is  the  relative  mass,  and  the  polarization  is 
equal  to  the  ionic  charge  times  the  relative  displacement  then  devided  by  the  unit  cell  volume. 
However,  in  the  perovskite  structure  there  are  three  different  types  of  ions,  hence,  a  more  realistic 
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model  would  be  a  three  body  system  model.  In  what  follows  we  will  use  a  one  dimensional  rigid 
ion  model  for  BaTiOs  as  an  example  to  illustrate  the  relationship  between  the  coefficients  in  Eu 
( 1 )  and  the  microscopic  quantities. 

According  to  the  structural  work  of  Shirane  et  al,*®  the  soft  optical  mode  in  BaTi03  consists 
of  the  relative  motion  of  Titanium,  Barium  and  Oxygen  (Fig.  la).  Because  the  center  of  mass  and 
the  center  of  charge  for  each  type  of  ions  coincide,  we  can  effectively  treat  this  system  as  a  three- 
body  system  in  the  lattice  dynamical  calculations.  For  convenience  the  ion  groups  are  labeled  as 
follows  (see  Fig.  lb):  Ba  —  1;  Ti  —  2;  and  3  O  — 3. 

The  potential  energy  represents  the  k  =  0  mode  for  this  three-body  system  is 

<I>  =  ^(Xi-X2)2-»-^(X2-X3)2  (19) 

In  order  to  derive  the  equations  of  motion  one  should  also  consider  the  Lorentz  field,  which  leads 
to  the  following  differential  equations: 


mi  Xia-Ki  (Xi-X2)  +  ^qiP 

(20  a) 

m2  X2  =  -  Ki  (X2  -  Xi  )  -  K2  (X2  -  X3  )  +^q2  P 

(20  b) 

m3  X3  =  •  K2  (X3  -  X2  )  +  ^q3  P 

(20  c) 

qi  +q2  +q3  =  0 

(21) 

P«(qi  Xi-fq2X2  +  q3X3)/a^ 

(22) 

where  ao  is  the  lattice  constant,  qi,  q2,  qs  are  the  charges  of  the  three  ion  groups,  P  is  the 
polarization. 

For  convenience  let  us  introduce  two  new  variables 
u  =  Xj  -X3;  V  =  X2  -  X3. 

Using  eqs.  (21)  and  (22),  and  the  new  variables  u  and  v  the  equations  of  motion  (2()a-c)  can  be 

simplified  to  become 

ii  s  a  u  -f  b  v  (23a) 

v  s  c  u  +  d  v  (23b) 
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Figure  1 .  illustration  of  the  ionic  displacement  in  BaTi03  from  ref. 
(b)  One  dimensional  modal  for  the  soft  mode 


where 

a=  -ii-* 
mi 


qi  f  SL  +  SL±Si ) 
3  £0  ^^3 


q2 


3£oa^ 


c  = 


m3  mi 

ki  ,  qi  (  q2  ,  qi  -Hq2  . 

"32  3eoa^  *"2  m3 


/  qi  ,  qi  +  q2 . 

^  m3  ' 


q2 


d  =  -(!^+ii)  + 

m2  m3'  3e„ai  ">2 


/  q2  ■  qi  +  q2  X 

^  ni3 


Eqs  (23a,b)  have  the  harmonic  oscillator  solution 
u  =  uo  exp(ja)t) 
v  =  voexp(j(Ot) 

and  the  eigenfrequency  to  is  given  by 


(24) 

(25) 


t4  =  ^  [  -a  -  d  ±  V (  a  +d)2  -  4  (ad  -  be)  ] 

In  any  given  mode  the  relative  displacement  u  and  v  are  proporuonal  to  each  other.  For  the 
coordinate  system  in  Fig.  1,  u  and  v  always  have  the  same  sign.  Assuming  one  of  the  modes, 
for  example  (O.,  becomes  soft  at  lower  temperatures  due  to  the  decrease  of  the  coupling  strength 
between  different  ions,  then,  their  magnitudes  have  the  following  relationship; 

vq  «  r\-  ( d  -  a  +  V(a  +  d)2-4(ad-bc) )  u©  (27) 

2  D 

From  Eqs.  (21)  and  (22)  the  polarization  P  is  given  by 

P  =  (qi  u  +  q2  V )/  ss(  qi  +  ;^  ( d  -  a  +  V(a  +  d)2-4(ad-bc) )]  u  /aj  (28) 

2  b 

Therefore,  in  this  mode 

Zes[qi+^(d-a  +  V(a  +  d  )2  - 4  ( ad  -  be ) )]  /aj  (29) 

2  D 

The  value  of  P  can  be  obtained  from  pyroelectric  measurements  and  u  may  be  calculated  horn 
X-ray  diffractions  of  the  low  and  high  temperature  phases,  hence,  Ze  can  also  be  obtained 
experimentally. 
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Since  v  and  P  are  linearly  proportional  to  u,  we  could  simplify  the  problem  by  constructing 
a  new  single  variable  potential  G=  (k/2)  which  gives  rise  to  the  following  equation  of  motion. 


=  ^  [  -(a+d)  -  V (a+d)^  -  4(ad-bc)  ]  u 


(30) 


where  k  has  the  dimension  of  force  constant  and  M  has  the  dimension  of  mass  according  to  the 
definition  of  a,  b,  c  and  d.  For  a  three  dimensional  system,  the  constructed  potential  which  leads 
to  Eq.(30)  will  have  the  same  form  as  Eq.  (1)  according  to  symmetry.  Because  u  is  also 
proportional  to  the  polarization  P  [eq.  (28)],  we  may  also  write  down  the  constructed  potential  in 
terms  of  the  polarization  vector  P,  which  becomes  the  Landau  potential  for  a  ferroelectric  system 
G=  (ce/2)  P2,  a  =  K/  (Z  e)2. 

For  longitudinal  vibrations,  we  have  to  add  the  depolarization  field  (-  P/Co)  to  the  equation 
of  motion  Eq.(20  a  <),  which  will  add  a  positive  contribution  to  the  eigenfrequency  preventing  it 
to  become  soft  like  the  transverse  mode.  Formally,  the  relationship  between  (Ol  and  cor  may  be 
written  as 

c^=o>f  +  ^  (30) 

where  A  is  a  positive  constant  reflecting  the  contribution  of  die  depolarization  field. 

IV-  Summary  and  wnclusidns 

It  is  shown  that  the  polarization  gradient  coefficients  in  the  Landau-Ginzburg  theory  can  be 
directly  related  to  the  dispersion  surface  of  the  soft  mode.  Therefore,  their  physical  meaning 
becomes  aj^iaient  The  correspondence  between  the  Landau  theory  and  the  lattice  dynamic 
potential  was  illustrated  duough  a  simple  one  dimensional  rigid  ion  model  for  BaTiOs  at  k=0. 
Since  the  polarization  and  the  relative  displacement  field  have  a  linear  relationship,  the  polarization 
gradient  coefficients  in  the  Landau-<}inzburg  theory  may  be  calculated  from  the  measured 
dispersion  curves  near  the  soft  mode.  For  cubic  symmetry  there  are  only  three  independent 
gradient  coefficients,  the  dispersion  anisotropy  of  the  soft  mode  can  be  determined  through 
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measuFements  along  the  three  principal  directions.  These  gradient  coefficients  in  principle  can  be 
obtained  through  inelastic  neutron  scattering  experiments.  However,  in  many  cases  these 
dispersion  curves  are  very  difficult  to  measure  because  of  the  high  transition  temperature.  To  my 
knowledge,  a  complete  set  of  the  dispersion  relations  do  not  exist  in  the  literatures  for  the  known 
ferroelectric  materials.  One  of  the  intentions  of  this  paper  is  to  re-emphasize  the  importance  to 
measure  these  dispersion  curves  which  can  be  used  for  the  study  of  domain  walls  in 
ferroelectrics. 

An  interesting  point  should  be  also  mentioned  is  the  possibility  of  obtaining  the  dispersion 
surface  of  the  soft  mode  through  direct  measurements  on  the  domain  wall  profiles^^  because  the 
gradient  coefficients  can  be  extracted  from  fitting  the  measured  polarization  profiles  to  the  soliton- 
like  solutions  of  nonlinear  nonlocal  continuous  medium  theor'.(^>*^)  This  could  be  very  useful  to 
study  the  dispersion  surface  of  the  over  damped  soft  mode,  such  as  in  BaTi03,  which  can  not  be 
obtained  through  inelastic  neutron  scattering. 
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Abstract 

Election  microscopy  studies  in  lanthanum  doped  lead  dtanate  reveals  the  evolution  of  a 
spadal  modulation  in  the  magnitude  of  the  spontaneous  pdarization  with  the  increase  of  the 
lanthanum  dopant  On  the  incorporation  of  ~  2S  atom  percent  lanthanum,  the  conventional  dOTiain 
structure  becomes  ill-defined,  and  tweed  microstiuctures  are  observed.  The  structural  information 
can  be  associated  witii  the  change  from  normal  ferroelectric  to  difhise  ferroelectric  phase  transition 
behavior.  DifTerent  fiom  twin  structures,  these  modulated  structures  represent  a  new  type  cf 
polarization  variation  existing  within  a  single  domain.  Further  understanding  of  die  observed 
spatial  variation  in  pdarization  requires  structural  analysis  at  the  atranic  scale.  Hologrtphic 
electron  microscopy  is  pn^iosed  as  a  potential  tod  to  study  various  polarization  gradients  in 
ferroelectric  materials.  Understanding  the  spatial  variations  in  polarization  is  essential  to  more 
fully  oooqnehend  the  extrinsic  contributions  to  the  elastonlielectric  properties  in  ferroeiectrics. 
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Introduction 


Mesoscopic  structures  within  ferroelectrics  and  related  materials  have  important 
consequences  with  respect  to  the  macroscopic  elasto-dielectric  properties.  These  structures  exist 
on  a  scale  of  a  few  tens  to  a  few  thousands  of  angstroms  and  include  defect  structures  within  the 
lattice  as  well  as  the  polarizadon  domain  structures  associated  widi  the  ferroic  j^ase  transition.  In 
general,  there  are  two  contributions  to  the  elasto-dielectric  pn^)erties:  the  intrinsic  contribution, 
which  is  related  to  the  ferroelectric  (antifenoelectric)  atomic  structure,  and  the  extrinsic 
contribution,  which  is  associated  with  domains  and  defects.  In  technologically  important 
materials,  such  as  Pb(Zr,'n)03,  the  extrinsic  factors  can  contribute  as  much  as  70%  to  the  total 
elasto-dielectric  response  (see  Figure  1).  Therefore,  it  is  necessary  to  develop  a  greater 
understanding  (tf  ail  the  possible  defect  and  polarization  mechanisms  which  could  contribute  to  the 
extrinsic  elasto-dielectric  properties.  However,  a  c(»iq>rehensive  theoretical  description  of  the 
extrinsic  ctmtributioiis  is  currently  not  in  place. 

The  most  common  mesoscopic  structures  associated  with  ferroelectric  and  related  materials 
are  domains  and  domain  walls.  Domains  form  at  the  phase  transition  and  relate  the  low 
tetr^erature  phase  to  the  high  temperature  prototype  phase  via  certain  symmetry  constraints.  In  the 
example  of  ‘‘normal”  first-  or  second-order  ferroelectric  transitions,  each  domain  is  an  area  of 
uniform  polarization,  and  the  boundary  which  divides  two  domains  (i.e.  a  twin  structure)  is  known 
as  the  domain  wall  Ihe  domain  wall  is  a  regitm  of  distorted  crystal  structure  in  which  there  exists 
a  spatial  transitkm  of  die  pohtrizarionfirom  one  orientation  state  to  another. 

There  are  two  main  types  of  twin  structures  Oie  type  is  a  twin  with  inversion  symmetry  of 
the  polarization  but  in  which  the  strain  is  die  same  in  bodi  variants.  The  second  type  is  a  twin  of 
two  variants  with  different  orientation  for  both  polarization  and  strain.  Ferroeiectric  twins  are 
typically  of  the  head  to  tail  configuraticm.  There  are  reports  of  other  domain  configurations,  such 
as  head-to-head  types,  but  these  have  not  been  Mtensivdy  investigated.^'*^ 

The  fine  structure  of  die  ferroelectric  domain  walls  depends  on  a  number  of  inter-related 
parameters:  including  the  symmetry,  temperature,  order  of  the  phase  transition,  spontaneous 
polarizaticm,  and  the  electrosttictive  and  elastic  cooqiliances.  A  number  of  analytical  descriptions 
now  exist  to  describe  the  spatial  variatkm  of  die  order  parameter  in  a  ferroelectric  donuun  waiL^<^7 
However,  sana  of  the  parameters  required  by  the  theory  cannot  be  easily  acquired  with  current 
experimental  techniques,  and  so  there  is  a  need  to  develc^  new  expoimental  methods  to  study 
polarization  variations  in  these  noaterials. 

Defects  and  dopants  are  known  to  have  a  strong  influence  on  the  elasiOKlieiectric  properties 
of  ferroelectrics  and  related  materials.  Theoretical  treatmentt  of  die  role  defects  near  structural 
phase  transitions  are  usually  restricted  to  defect  densities  that  are  much  less  than  the  reciprocal 
correlation  volume  (- 10^^  cm*^).^  In  the  perovsldte  ferroelectrics  of  conomercial  interest,  such  a 
situation  is  almost  never  realized. 
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This  article  outlines  some  of  the  results  observed  by  diffraction  contrast  electron 
microscopy  in  perovskite-based  ferroelectrics.^0  From  these  results,  and  requirements  of  the 
theoretical  developments,  there  is  an  urgent  need  to  experimentally  investigate  the  polari/ation 
gradients,  both  within  the  domain  region  and  in  the  region  of  the  domain  wall  Electron 
holography  is  discussed  as  a  technique  potentially  capable  of  quantifying  the  polarization  gradients 
in  these  materials. 

Results  and  Discussion 


Doping  in  Pb(Zr.Ti)03-based  materials  by  lanthanum  is  used  as  a  means  to  soften  the 
switching  characteristics  of  piezoelectric  materials.^!  Additionally,  the  incocpoiation  of  lanthanum 
facilitates  the  fabrication  of  transparent  ceramics  for  optoelectrtmic  applicatitms.^^  In  general, 
doping  with  lanrhanum  has  a  significant  influence  on  many  of  the  elasto-dielectric  prq)etties.  For 
sufficiently  high  levels  of  doping  in  Zr-rich  PZT  compositions,  this  leads  to  the  observation  of 
diffuse  phase  transition  behavior  having  strcmg  dielectric  dispersion.  Fenoelectrics  with  this 
behavior  are  generally  referred  to  as  relaxors,  and  are  of  technological  inqxxtance  owing  to  their 
unique  electrostrictive,  ciqxicitive,  and  optoelectronic  propenies.  The  domain  structures  of  relaxor 
(PbJji)(Zr,Ti)03  (PLZT)  ceramics  are  difficult  to  study  using  transmission  electron  microscq)y. 
However,  by  carefiiUy  coding  a  8  JZ^O/30  coa9ositk)n.  a  nticrodonaain  contrast  could  be  detecmd, 
as  shown  in  Hgure  2(a).  Under  the  irradiariOT  of  die  electron  beam,  die  domain  structure  is 
unstable.  By  agitating  the  structure  through  focusing/defocusing  the  beam,  the  domain 
ctmfiguration  transforms  to  a  more  stable  ordered  structure.  Figure  2(b).  It  is  believed  that 
dietmally-induoed  stresses  switch  the  mictodomain  structure  to  a  new  domain  configuration. 

The  end-member  of  die  PZT  soUd  sdutkm,  PbHOs,  has  the  hi^iest  tcanridon  tenqierature 
(Tc  »  490  and  the  largest  strain  [(c/a  -1)  «  6.S  %]  widiin  the  petovskite  family.  This  makes 
PbTiOs  an  ideal  material  to  stwfy  by  transmission  election  ndcroscopy.  Doping  PbTiOs  with 
lanthannm  (WT)  wdnces  both  the  phase  transition  temperature  and  the  characteristic  discontinuity 
of  the  first-order  transition.  A  systematic  study  of  the  structural  effect  of  landianum  on  the  polar 
Hntnain  structuTB  in  PLT  ccruiiics  reveals  die  development  of  a  strain  texture  within  the  normal 
dmnains,  Hgure  3(a),(b),(c).  Using  diffiraction  contrast  invisibility  conditions,  we  were  able  to 
deduce  that  the  texture  is  the  result  of  a  non-uniform  ^lontaneoas  deformation  along  the  c-axis 
within  die  /inmain,  As  shown  in  Hgure  3(a),  for  a  sample  doped  widi  1  atom  %  La,  there  is  no 
evidence  of  a  texture.  At  lanthanmn  concentration  is  increased  from  between  5  and  10  atom  %. 
texturing  appears,  and  this  becomes  progressively  more  pronounced  with  increasing  dopant 
concentration.  When  the  dopant  concentration  readies  2S  atom  %,  a  noniud  domain  structure  is 
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not  identified,  and  a  full  cross-hatched  or  "tweed"  domain  strucmre  is  observed  below  che 
transition  temperature,  Hg.  3d.  Similar  structures  have  been  observed  in  ferroelastic  systems  such 
as  YBa2(Cuo.9Feo.03)307-5  and  Mg-  Coniicrite.i^**5  Inhomogeneous  polarization  distributions 
are  not  unique  to  the  PLZT  system,  but  also  exist  in  the  complex  lead  Pb<B3'0O3  perovskite 
systems.  In  these  systems,  the  intermediate  scale  B-site  cation  <»dering  is  the  source  of  the 
polarization  modulation.  In  order  to  further  our  understanding  of  the  polarization  variation  in 
these  materials,  we  need  to  develop  a  technique  to  quantify  the  polarizadon  gradients  and  defect 
structures.  In  this  regard,  die  potential  of  the  electron  holography  technique  is  discussed  below. 

Elccttnn  Holognphy 

The  idea  of  using  coherent  electrons  in  electnm  microscopy  was  proposed  in  1949  by 
Gabor  in  an  attempt  to  extend  the  limits  of  electron  microscope  resolution.  However,  the 
realization  of  electron  holognqihy  was  achieved  cmly  in  the  1980*s  owing  to  the  development  of  a 
coherent  field-enusskm  electron  beam.  Commercial  instruments  for  electron  holography  have  be«i 
developed  by  Ifitachi  Lid.  and  Philips  but  have  only  recently  become  available.  The  principle  of 
electron  holography  is  similar  to  that  of  optical  holography,  in  which  the  die  phase  and  amplitude 
of  die  electrao  beam  are  lecoided  simultaneously.  The  addition  of  phase  shift  infixmation  which  is 
highly  sensitive  to  local  changes  makes  electrra  htdogr^il^  a  more  attractive  method  compared  to 
conventional  election  microscopy  tediniques.  There  have  been  a  variety  of  applications  for  this 
new  technique  starting  since  1980,  especially  in  the  study  of  magnetic  domains  and  fluxons  in 
superconducting  materials, 

Recently  the  possibility  of  using  die  holognqdiy  technique  to  study  ferroelectric  domain 
walls  and  odier  defect  structures  in  fenoelectiics  was  recognized.  Some  encouraging  results  have 
bera  reported  on  the  profiles  of  domain  walls,  as  shown  in  Figure  The  kink-like  electron 
interference  fringe  pattern  doady  resembles  the  qiace  profile  of  the  polarization  vector  across  a 
domain  wall  as  predicted  by  the  Landau^jinzborg  model.^>d  Although  a  conqilete  theoretical 
description  of  the  fringe  profile  in  Ref.  18  is  not  cunently  available,  the  fhct  that  dm  electric  field 
variation  can  be  probed  on  a  scale  lessthanlA  is  both  exciting  and  promising. 

It  has  been  demonstrated  that  the  electron  holography  technique  may  also  be  used  to  study 
the  location  of  aliovalent  deviants  inside  crystal  structures  through  the  perturbed  local  electric 
fidds.^  As  shown  in  Fig.  5a  the  fringe  bifurcations  oocuned  across  die  domain  walL  Potential 
ctmtours  reveals  diere  ate  charge  centers  attracted  to  the  domain  wall  (see  Figure  5b).^  This 
information  may  lead  to  a  significant  advance  in  the  understanding  of  the  effect  of  dopants,  and 
may  shed  new  tidit  on  the  study  of  interactions  between  the  dopants  and  domain  walls  and  other 
pdarization  modulttioos  as  described  above. 
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Quantitative  study  of  the  polarization  profiles  can  have  a  significant  impact  on  the 
fundamental  understanding  of  fenoelectrics.  Once  the  relationship  between  the  fringe  variation  and 
the  polarization  space  profile  is  established,  one  can  obtain  the  polarization  gradient  coefficients 
through  back  fitting  the  observed  domain  profiles  to  the  theoretical  results  on  dmnain  walls  J  These 
gradient  coefficients  are  a  measure  of  the  nonlocal  coupling  strengtii.  Using  lattice  dynamics,  one 
may  correlate  these  gradient  coefficients  to  the  dispersion  surface  near  the  soft  mode  of  the 
paraelectric-ferroelectric  phase  transition.^!  Hence,  the  electron  holography  technique,  together 
with  the  continuum  model  described  in  Ref.  6,  can  potentially  provide  a  methodology  to  study  the 
characteristics  of  the  over-danced  soft  mode  in  systems  such  as  BaTiQs,  which  could  not  be 
directly  probed  through  inelastic  neum»  scattering. 

As  a  new  technique,  many  problems  still  exist  in  tire  electron  holography,  especially  with 
regard  to  the  interpretation  of  the  observed  fringes.  In  principle,  the  total  phase  shift  of  high 
energy  coherent  electrons  passing  throu^  a  friroelecttic  thin  sancleilPy  be  calculated  frmn  the 
following  equatioi^ 

♦(xo.yo)*j^|v(xo,yo.*)dz,  (1) 

■f- 

wlrere  X  is  the  electron  wavelengdi,  Xq  and  yo  define  die  point  on  die  thb  sample,  £is  the  electron 
energy,  and  V(xo,y<hZ)  represents  die  electrical  potential  experienced  by  the  traveling  electrons. 
However,  V(xo,yo,z)  represents  the  total  potential,  and  it  is  quite  difficult  to  delineate  contributions 
from  the  **bouiur  charges  (relevant  to  the  polarizatkm)  and  the  “free**  charges  (relevant  to 
condensation).  We  believe  diis  is  the  main  reason  for  die  inconsistencies  encountered  in  the 
current  studies  fenoelectiics  using  electron  hoiognqihy.^  More  theoretical  analysis  of  the 
interptetadon  (tf  the  htdogrqdty  results  is  in  order. 

Conclusions 

Observiliaot  by  convencioaal  transinissioo  etoctron  tmcroscopy  techniques  on  ferroelectric 
and  related  materials  reveal  a  varies  of  polarization  modulations  uhich  can  be  induced  when  there 
exists  coupling  of  die  primary  order  parameter  to  symmetry  breaking  defects.  I^om  the  evolution 
of  the  modulated  structures  and  domain  structures,  one  can  see  some  link  between  these 
mesoscopic  siruenues  and  the  extrinsic  elasto-dielectric  properties. 

A  new  electron-microscqpy  technique  using  coherent  electrons  known  as  electron 
holography  opens  up  opportunities  in  die  study  ci  domain  walls  and  defect  structures.  In  this 
technique,  phase  shifts  can  be  correlated  to  local  variations  of  die  electrottatic  fields  within 
materials.  Further  development  of  this  new  technique  in  tiie  study  of  fetrodectrics  may  help  us  to 
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gain  new  insight  into  the  mechanisms  of  extrinsic  contributions  to  the  macroscopic  elasto-dielectric 
properties. 
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Figures 


Temperature  dependance  qf  the  dielectric  permitdvity  of  doped-PZT  ceramics.  The 
theoiedcal  value  from  the  Landau-Devonshiie  theory  is  also  shown. 

(a)  Bright  field  image  of  liquid  nitrogen  cooled  PLZT  S.2/70/30  relaxor  ferroelectric 
revealing  a  microdomain  structure.  ( b)  An  in  situ  switched  pseudo-domain 
structure  of  (a). 

Bri^t  field  image  of  domain  structures  in  a  sdid  solution  series 
(Pbl-ax/zLaOTlOs  with  four  different  La  contents. 

An  electron  hologram  a  90**  domain  wall  in  BaTiOj,  the  fringe  bending  is  related 
to  the  polarization  difference  across  the  domain  wall  (ref.  19). 

(a)  An  electron  hologram  showing  anomalous  fringe  bifurcations  in  a  90**  domain 
wall  in  BaTiOs.  (b)  Electron  interfetogram  of  the  same  area  o{  (a)  shows  charge 
defea  centers  within  the  wall  (courmsy  of  Drs.  D.  Joy  and  X.  Zhang).^ 
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MICROSTRUCTURAL  CHARACTERISTICS  AND  DIFFUSE  PHASE 
TRANSITION  BEHAVIOR  OF  LANTHANUM-MODIFIED  LEAD  TITANATE 
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Abstract  The  introduction  of  structural  disorder  into  perovskite  ferroelectrics  leads 
to  the  observation  of  diffuse  phase  transition  behavior.  In  this  paper,  the 
microstructural  characteristics  and  phase  transition  behavior  of  lanthanum- 
modified  lead  titanate  are  discussed.  It  is  shown  that  the  diffuse  nature  of  the 
transition  is  connected  witti  the  appearance  of  ferroelectric  domain  structures 
exhibiting  texture  on  d\e  mesoscopic  ( «  10  nm)  length  scale.  The  theoretical 
implications  of  these  results  are  briefly  discussed. 

INTRODUCTION 

Diffuse  phase  transitions  in  perovskite  ferroelectrics  occur  as  a  consequence  of 
some  level  of  structural  disorder  that  breaks  the  translational  invariance  of  the 
lattice.  The  type  of  diffuse  transition  observed  depends  both  on  the  nature  and  scale 
of  the  structural  disorder.  Impurities,  point  defects,  extended  defects,  incomplete  or 
inhomogeneous  cation  ordering,  macroscopic  fluctuations  in  chemical  composition, 
core/shell  structures  etc.  can  all  lead  to  smearing  of  the  ferroelectric  phase 
transition.  The  physical  origin  of  the  smearing  and  the  behavior  of  the  ferroelectric 
properties  in  each  particular  case,  however,  may  be  very  different.  In  chemically 
complex  compounds  and  solid  solutions,  several  types  of  structural  disorder  are 
often  present  simultaneously.  To  establish  structure-property  relations  in  these 
materials  therefore  requires  characterization  at  all  the  relevant  length  scales, 
including  the  macroscopic  (crystal  symmetry,  chemical  homogeneity,  core/shell 
structures),  the  microscopic  (ferroelectric  domain  structures,  extended  defects), 
the  mesoscopic  (cation  order  domains,  texture  within  ferroelectric  domains  ),  and 
the  atomic  (point  defect  chemistry,  local  environment). 

For  ferroelectrics  such  as  lanthanum-modified  lead  zirconate-titanate  (PLZT), 
the  structural  disorder  arises  as  a  consequence  of  the  compiex  defect  chemistry 
associated  with  the  aliovaient  A-site  substitution  of  lanthanum^.  Keizer,  Lansink 
and  Burggraaf^  have  shown  that  the  nature  of  the  ferroelectric  transition  In  ceramic 
specimens  of  the  end  member  Pbi-ayLayTiOs  (PLT)  changes  smoothly  from  that 

expected  in  a  sharp  first-order  phase  change  to  that  of  a  diffuse  transition  as 
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the  lanthanum  concentration  is  increased.  Figure  1  shows  that,  when  characterized 
by  the  exponent  yin  the  generalized  Curie-Weiss  law,  the  anomalous  dielectric 
behavior  onsets  at  y  «  0.05  and  becomes  typical  of  diffuse  transitions  (/.  e.  y  »  2) 
when  y  >  0.23.  Although  the  numerical  values  of  y  depend  somewhat  on  the  lead 
elimination  factor  (a)  as  well  as  on  the  grain  size,  the  general  trend  towards 
increasingly  diffuse  behavior  with  lanthanum  substitution  is  intrinsic. 

In  order  to  better  understand  how  structural  disorder  evolves  and  relates  to  the 
dielectric  properties  in  PLT,  we  have  recently  undertaken  transmission  electron 
microscopy  (TEM)  and  powder  x-ray/neutron  studies,  the  details  of  which  are 
reported  elsewhere^'S.  in  this  paper,  we  draw  on  these  results  and  speculate  on 
some  possible  interconnections  between  the  microstructure,  dielectric  properties, 
and  the  phase  transition  behavior. 


High  resolution  X-ray  diffraction  studies^  of  sol-gel  derived  powders  have  shown 
that  the  addition  of  lanthanum  to  lead  titanate  results  in  broadened  diffraction 
profiles  showing  a  marked  profile  asymmetry.  The  degree  of  asymmetry  can  be 
quantified  by  fitting  ttte  profiles  to  split  Pearson  distributions^.  The  split  Pearson 
function  gives  two  values  for  the  full  width  at  half  maximum,  p*-  arnf  5^,  which  are 
characteristic  of  the  peak  shape  on  the  low  and  high  angle  sides  of  the  peak  maxknum, 
respectively.  Figure  2  shows  that  both  ttte  peak  asymmetry  and  the  peak  breadth 
renormalize  as  a  function  of  lanthanum  concentration.  The  asymmetry  (5^/5^) 
changes  markedly  between  y  >  0  and  y  «  0.01,  but  then  decreases  and  becomes 
constartt  for  samples  with  y  >  0.05.  At  the  same  time,  the  profile  breadth  decreases. 
Corresponding  to  the  changes  in  the  diffraction  profiies,  the  crystai  tetragonaiity 
(c/a)  begins  to  deviate  from  the  iinear  composition  dependence  as  expected  based  on 
Vegard's  iaw.  These  observations  cannot  be  easiiy  expiained  by  macroscopic 
variations  in  chemicai  composition,  nor  by  particie  size  effects. 

Aitemativeiy,  recent  TEM  studies^  (Figure  3)  reveai  that  the  deveiopment  of 
sub-domain  texture  in  PLT  materiais  cioseiy  paraifeis  the  changes  observed  both  in 
the  X-ray  studies  and  in  the  dielectric  behavior.  For  samples  with  y  <  0.05,  no  sub- 
domain  texture  is  observed,  and  the  transition  behavior  is  sharp.  For  samples  with 
y  z  0.05,  sub-domain  texture  becomes  apparent,  and  increases  in  degree  with 
lanthanum  concentration.  As  the  sub-domain  texture  becomes  more  prorraunced,  the 
dielectric  behavior  near  the  transition  becomes  correspondingly  more 
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Figure  3.  'Ranmissloii  electron  mtcrographe  of  Pbi^^^^TlOs: 
W  y  ■  0.01  (fa)  y  ■  0.05  (cj  y  ■  0.10  and  (<9  y  ■  oJs. 
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anomaiotw.  For  compositions  with  y  >  0.25.  the  microstructure  degenerates  into  a 
cross-hatched  or  tweed  texture.  For  these  compositions,  the  dielectric  behavior  is 
typical  of  that  for  diffuse  phase  transitions.  Note  that  these  textures  disappear  on 
heating  above  the  transition  temperature,  and  so  cannot  be  identified  with  exsolution 
lamellae. 

The  tweed  structures  observed  for  the  higher  lanthanum  concentrations  would 
appear  to  be  precursors  for  the  mottled  and  poorly  defined  vestiges  of  sub-domain 
texture  as  observed  for  PLZT^.  The  Zr-rich  PLZT  compositions  exhibit  strongly 
frequency-dispersive  diffuse  transitions  similar  to  those  of  the  complex  B-site 
perovskites  such  as  PbMgi/3Nb2/303  (PMN)  and  related  compourKis.  In  the  later 
case,  however,  it  is  the  cation  order  domains  persisting  on  a  -10  nm  length  scale 
that  lead  to  the  breakdown  of  the  ferroelectric  domain  structure^'^^.  The  similarity 
in  mesoscopic  domain  texture  of  the  PLZT  and  PMN  type  materials^ ''*^3  is 
intriguing,  since  at  the  atomic  level,  the  origin  of  the  structural  disorder  in  these 
two  materials  is  very  different. 

With  regard  to  the  nature  of  the  phase  transition  in  PLT,  the  first-order 
character  also  follows  the  changes  in  structural  properties  as  revealed  in  the  X-ray 
data.  The  development  of  a  Landau-Oevonshire^^  formulation  for  the  phase 
transition  in  PLT  requires  the  expansion  coefficients  to  depend  on  spatial  variables, 
which  results  in  a  distribution  of  transition  temperature.  As  a  macroscopic  average, 
one  expects  that  the  first-order  discontinuity  of  the  transition  will  be  drastically 
reduced  as  the  transition  becomes  diffused.  The  Landau  theory  is  valid,  however, 
only  when  the  sub-domain  modulations  are  weak,  otherwise  the  gradient  related 
(Ginaburg)  terms  must  be  indudad.  Uskig  the  methodology  described  previousiy^S, 
the  coefficients  in  the  Landau-Oevonshire  expansion  were  estimated  for  compositions 
with  y  <  0.05.  where  judging  from  Figure  3.  this  condition  is  approximately 
satisfied.  As  shown  in  Figure  4,  the  coefficient  of  the  quartic  term  in  the  elastic 
Gibbs  function  (an^)  increases  rapidly  toward  zero  between  y  «  0  and  y  ■  0.01  but 
then  changes  at  a  much  lower  rate.  Utis  increase  reflects  the  apparerrt  loss  in  the 
first-order  character  of  the  transition  as  the  structure  relaxes  to  the  defects  and  the 
long  range  order  is  disrupted.  The  coefficient  of  the  sixth-rank  term  also 
renormaHzes,  but  to  a  much  lesser  extent.  Large  changes  in  the  coefficient  of  the 
quartic  term  were  also  observed  for  PLZT  compositions^^. 

For  the  compositions  with  y  >  0.05,  the  Ginsburg  terms  must  be  included  in 
the  formaltom.  The  inclusion  of  the  gradient  energy  can  describe  both  the  domains 
and  thair  modulations.  Periodic  domain  structures  can  be  excited  for  a  finite 
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system  du«  to  the  competition  between  noniinearity  and  nonlocaiity.  Moduiated 
structures  and  tweed  textures  can  be  described  in  terms  of  crest  riding 
periodons^^'^^.  However,  the  moduiated  structures  such  as  those  in  Figure  3c  do  not 
seem  to  follow  the  expected  temperature  variation  pattern  of  the  crest  riding 
periodon  excitations,  which  are  unstable  and  can  survive  only  in  a  finite 
temperature  range  near  the  transition.  One  possible  explanation  could  involve  an 
interaction  of  the  defects  with  the  periodons.  These  excitations  might  be  stabilized 
by  the  lanthanum/vacanciea  if  the  defects  are  appropriately  aligned.  The  alignment 
can  be  driven,  on  the  other  hand,  by  the  field  gradient  created  at  the  walls  of  the 
periodons.  Consistent  with  this  notion,  the  modulations  have  preferred  orientations 
(although  not  perfect),  and  are  nearly  periodic,  as  shown  in  Figure  3.  Further 
studies  will  be  required  to  confirm  this  picture. 

At  the  higher  lanthanum  concentrations  (y  >  0.25)  the  microstructure  breaks 
down  on  still  finer  scales  and  more  degenerate  states  are  created.  This  will 
effectively  increase  the  dielectric  response  of  the  system,  as  observed.  For  these 
materials,  however,  arguments  based  on  the  continuum  theory  as  summarized  above 
may  not  be  used  because  the  density  representation  (both  in  terms  of  energy  and 
polarizalion)  is  no  longer  valid.  Instead,  dtecrete  models  must  be  developed.  The 
results  described  above  can  be  represented  schematicaliy  as  shown  in  Figure  5, 
where  the  observed  relationships  between  the  first-order  character  of  the 
transition,  the  crystal  tetragonality,  the  sub-domain  texture,  and  the  dielectric 
properties  are  shown.  Here  it  is  worth  mentioning  that,  in  the  Landau-Oevonshire 
formalism,  close  irrter-relationships  exist  between  the  curvature  of  the  energy 
surface,  the  order  of  the  transition,  and  the  symmetry  of  the  ferroelectric  phase.  In 
this  corwiection,  it  is  of  interest  to  consider  how  the  scenario  depicted  in  Figure  5 
may  change  as  zirconium,  which  lowers  the  tetragonality  of  the  system  and  further 
flattens  the  energy  surface.  Is  added  to  the  system.  Based  on  the  above  discussion,  we 
might  expect  that  lanthamim  will  become  more  effective  at  breaking  down  ttte  domain 
structure  as  the  rhombohedral  side  of  the  PLZT  diagram  is  approached,  leading  to  a 
large  number  of  degenerate  states  and  additional  contributions  to  the  dielectric 
response  even  at  modest  lanthanum  concentrations. 

SUMMARY 

The  breakdown  of  conventional  ferroelectric  phase  transition  behavior  in 
ianthmum-modified  lead  titanate  has  been  correlatsd  wW)  the  appearance  of 
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domain  structuras  exhibiting  texture  on  a  nanometer  scaie.  It  was  suggested  that 
this  type  of  mesoscopic  sub-domain  texture  is  a  common  feature  of  ferroelectrics 
showing  diffuse/frequency  dispersive  transitions.  The  theoretical  implications  of 
the  relationships  existirtg  between  the  curvature  of  the  energy  surface,  the  order  of 
the  trarwition,  the  symmetry  of  the  ferroelectric  phase,  the  development  of  the  sub- 
domain  texture,  arMI  the  resulting  dielectric  properties  were  discussed. 
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I.  INTRODUCTION 

This  is  the  progress  report  for  the  first  year  of  a  two  year  program  sponsored  by  the  Air 
Force  Office  of  Scientific  Research  on  "Domain  Processes  in  Ferroelectric  Ceramic".  Domain 
processes  contribute  a  large  poition  in  the  dielectric  and  piezoelectric  activities  in  ferroelectric 
ceramics,  a  better  understanding  of  this  subject  can  be  of  great  help  in  future  materials  development 
tor  acoustic  devices  and  smart  structures.  The  matoiai  studied  in  the  first  year  investigation  is  the 
PbZiOs-PbTiOs  solid  solution,  the  so  called  PZT  ferroelectric  ceramic,  which  has  been  the  primary 
transduce  material  fOT  the  past  thirty  years.  Although  PZT  is  a  well  known  transducer  material  and 
has  been  used  widely,  there  has  not  been  any  systematic  theoretical  study  in  die  past,  which 
unambiguously  characterizes  this  system,  especially  lacking  the  knowledge  on  the  mechanism  for 
the  strong  piezoelectric  activities.  Only  in  the  last  few  yean  people  have  starmd  to  realize  that  the 
major  conoibutor  in  the  piezoelectric  effect  is  fiom  the  domain  processes  and  have  made  a  start  to 
distinguish  and  separate  the  intrinsic  and  extrinsic  contributions.  Our  focus  in  this  project  is  to 
undentand,  firom  a  theoretical  point  cf  view,  the  extrinsic  contributions.  In  the  first  year  of  this 
Air  Force  sptxisored  program,  we  have  carried  out  a  critical  study  on  the  PZT  system  and  found 
many  ill-defined  ctmcepts  and  definitions.  Anoong  these  are  the  concept  of  the  morphotropic  phase 
boundary  (MPB),  die  distribution  of  coexisting  phases  in  a  complete  solid  solution  system,  the 
rii-«uinctMMi  between  domain  switching  and  field  induce  phase  transition,  and  the  relationship 
between  ynfipiMi«e  boundaries  and  ferroelectric  domain  walls  in  die  rhombohedral  II  phase.  These 
problems  have  been  addressed  in  die  first  year  and  the  details  ate  given  in  die  following  sections 
n,  in  and  IV,  respectively.  A  recently  funded  University  Initiative  Research  (URI)  program  by  the 
Office  of  Naval  Research  provided  a  solid  back  up  on  the  experimental  part  oS  diis  study. 
Collaborative  work  between  the  two  programs  is  already  in  progress  (see  section  IV). 
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n.  STUDY  ON  THE  MORPHOTROPIC  PHASE  BOUNDARY  IN  PZT 

The  most  useful  PZT  compositions  are  those  near  the  so  called  Morphotropic  Phase 
Boundary  (MPB).  The  MPB  is  defined  as  a  compositional  phase  boundary  at  which  the  free 
energies  for  the  rhombohedral  and  tetragonal  phases  are  equal.  It  is  well  known  that  many 
physical  properties  peak  at  the  MPB,  but  the  reason  for  this  is  not  well  understood.  In  addition, 
the  fitee  energy  is  not  a  directly  measurable  quandty,  in  order  to  study  the  MPB  we  need  to 
conrelate  the  equal  fiee  energy  concept  to  some  measurable  quantities.  In  the  past,  thennodynamic 
concepts  devel<^)ed  for  liquid-solid  transition  and  liquid  (or  gas)  mixtures  were  borrowed  to 
address  the  MPB  and  the  distribudon  of  structural  phase  mixture  in  the  PZT  system.  As  a  result, 
some  inconsistencies  and  confusions  are  created.  In  light  of  this  situadon,  we  have  conducted  a 
systematic  theoretical  study  on  the  MPB  in  the  PZT  system  during  the  first  year  of  this  project 
A.  The  Definition  of  MPB 

As  panted  out  in  Appendix  I  that  die  MPB  is  conventionally  represented  by  a  neariy 
vertical  line  on  the  phase  diagram  (see  figure  1  in  Appendix  I).  In  the  classical  book  on 
piezoelectric  ceramic  by  Jaffe  Cook  and  Jaffe,(l)  the  MPB  line  is  measured  by  X-Ray  diffraction 
technique.  Thenearly  vertical  line  on  the  phase  diagram  in  the  vicinity  of  52/48  Zi/11  ratio  is 
obtained  acooiding  to  the  criterion  ci  equal  volume  fraction  of  die  tetragonal  and  rhombohedral 
phases.  In  a  liquid  (or  gas)  mixtuie,  the  equal  volume  criterion  can  be  derived  from  thermodynamic 
principles  for  an  equilibrium  system.  However,  for  the  structural  phase  mixtuie  (which  may  not 
evo)  be  in  equilibrium),  the  two  phases  derived  from  a  common  parent  phase  have 

different  degree  of  freedom  in  the  configurational  space.  It  is  shown  in  Appendix  I  dtat  the  MPB 
defined  according  to  equal  free  energy  criterion  should  not  have  equal  volume  fraction  for  the  two 
coexisting  phases  if  the  geometrical  constraints  are  imposed.  Because  the  two  low  temperature 
phases  are  derived  from  the  same  cubic  phase,  in  order  to  calculate  the  volume  fractions  we  must 
stutfy  the  actual  transition  process  and  include  the  structural  constraints  in  the  transition  probability 
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calculation.  It  is  shown  (see  Appendix  I)  that  the  volume  ratio  for  the  riicanbohedral/tetragonal 
phases  should  be  60:40  at  the  MPB  instead  50:50  used  before. 

B.  Phase  Mixing  in  Complete  Solid  Solution  System 

As  shown  in  the  phase  diagram  (see  figure  1  of  Appendix  I)  PZT  is  a  complete  solid 
solution  system  in  the  cubic  phase.  There  is  a  paraelectric-ferroclectiic  phase  transition  which  occur 
at  around  370  ^  C  for  the  MPB  composition.  In  such  low  temperature  regime,  there  is  no  chemical 
driving  force  for  phase  separation  to  occur  in  this  system.  At  the  transition,  the  system  only 
experiences  a  displacive  structural  distortion  which  produces  a  spontaneous  polarization  due  to  the 
relative  shift  of  the  positive  and  the  negative  charge  centers.  X-ray  diffraction  shows  that  the 
symmeay  of  the  low  temperature  ferroelectric  phase  could  be  either  tetragonal  or  rhombohedral, 
depending  on  the  Zi :  Ti  ratio  in  the  system.  Near  the  MPB  composition,  the  low  temperature  phase 
is  a  mixture  of  the  tetragonal  and  rhombohedral  phases.  In  an  attempt,  die  fraction  of  each  low 
temperature  phase  was  described  by  the  lever  rule  and  the  two  edge  compositions  were  fitted  from 
the  experimental  data.C2*3)  Although  the  fitting  appears  to  be  quite  good,  the  edge  compositions, 
which  specify  the  width  of  the  coexistence  region,  can  not  be  uniquely  determined,  width  of  the 
coexisting  regimi  ranging  from  0  to  20  mole  percent.  It  is  found  that  in  general  ceramics  of  smaller 
grain  size  have  larger  coexistence  region  and  ceramics  of  larger  grain  size  have  smaller  coexistence 
region. 

We  have  looked  into  diis  problem  in  more  detail  and  found  that  the  description  the  phase 
coexistence  in  FZT  using  the  lever  rule  is  actually  in  ctmtradiction  widi  the  nature  of  the  complete 
solid  solutioo  syssem,  because  the  two  edge  compositions  in  die  lever  rule  actually  specify  a 
solubility  gq).  Since  such  solubility  gap  does  not  exist  in  the  PZT  system,  naturally  one  can  not 
define  these  two  edge  compositions  using  the  lever  rule.  Considering  die  fact  that  the  pataelectric- 
fertoelectric  phase  transition  is  second  order  for  PZT  at  conqiositions  near  the  MPB,  we  propose 
diat  the  phase  coexistence  in  this  system  is  a  result  of  a  "frozen  in"  metastable  phase  from  thermal 
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fluctuations  at  the  paiaelectric-fenoelectric  phase  transition  (see  Appendix  II).  A  method  is 
introduced  to  calculate  the  fractions  of  the  two  "frozen  in"  phases,  which  takes  into  account  both 
the  structural  constraints  and  the  statistics.  The  theoretical  calculations  agree  well  with  the 
experimental  results.  More  importantly,  the  new  model  does  not  contain  the  two  fixed  edge 
compositions  as  in  the  lever  rule,  instead,  two  new  variables  are  introduced:  one  is  the  MPB 
composition  which  is  a  fixed  quantity,  and  another  is  the  width  cf  the  coexistence  region  which  is 
inversely  proportional  to  the  volume  of  the  particles  (or  grains).  The  new  model  not  only  is 
compatible  with  the  nature  of  solid  solution  system  but  also  can  explain  the  grain  size  dependence 
of  the  coexistence  region  observed  experimentally.  Although  the  quantitative  test  of  the  new  model 
awaits  more  refined  experimental  input,  we  believe  that  the  conceptual  breakthrough  is  significant 
and  could  trigger  further  development  on  this  subject 

m  STUDY  THE  PHASE  TRANSmON  IN  RHOMBOHEDRAL  PHASE  P2T  AND 
ANTIPHASE  STRUCTURE  IN  R3c  PHASE 

As  part  of  the  first  year  effort  we  have  initiated  a  joined  investigation  on  the  PZT  system 
with  the  Navy  sui^xnted  URI  program .  The  focus  is  on  the  R3m R3c  phase  transition  in  the 
rhombohedtal  ferroelectric  phase  of  PZT.  This  phase  transition  is  not  well  understood  in  terms  of 
the  Landau-Devonshire  formalism.  In  addition  to  the  soft  zone  center  mode,  there  exists  anotiier 
soft  zone  boundary  mode  which  involves  the  tilting  of  the  oxygen  octahedra  in  the  perovsldte 
structure  and  causes  unit  cell  doubling.  We  have  made  a  detailed  structural  study  using 
TraiisiitissioiiBiectiOfi  Microscopy.  One  of  the  most  interesting  findings  is  that  many  antqrhase 
boundaries  ate  generated  in  the  R3c  phase  and  these  antiphase  boundaries  are,  in  most  cases^  not 
confined  by  the  crystal  symmetry  (see  Appendix  III).  The  antiphase  boundaries  with  their 
orientation  compatible  with  the  crystal  symmetry  have  been  addressed  by  Cao  and  Barsch,(^  but 
the  observed  antiphase  boundaries  in  the  PZT  system  seem  to  be  not  confined  by  syimnetry  at  all. 
We  suspect  that  there  ate  two  different  kinds  of  antiphase  boundaries  in  the  R3c  phase  of  PTZ,  in 
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one  kind  the  boundaries  are  orientated  along  the  polarization  direction  and  the  other  kind  with  their 
orientation  perpetKlicular  to  the  polarization.  The  observed  'uitiphase  boundary  may  be  the 
combination  of  the  two.  Modeling  of  the  antiphase  boundary  along  this  line  is  currendy  in 
progress. 

PZTs  of  the  rhombohedral  symmetry  is  one  of  the  best  pyroelectric  materials. 
Unfortunately,  the  existence  of  this  R3m  R3c  phase  transition  produces  a  second  kink  in  the 
polarizaticm  versus  temperature  plot,  which  effectively  limits  the  temperature  range  in  real 
applications.  The  understanding  on  this  transition  process  and  the  associated  domain  structures 
could  help  to  overcome  this  hindrance. 
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VI.  PROGRAMS  IN  PROGRESS 

Several  important  issues  have  attracted  our  attention  through  the  first  year  investiganon. 
However,  due  to  the  limitation  of  manpower,  it  is  not  possible  to  address  all  of  them.  Currently, 
based  on  the  already  achieved  progress,  we  are  pursuing  the  following  tc^ics: 
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a)  Study  the  periodic  structure  of  domains  in  ferroelectric  systems.  This  includes  theoretical  study 
on  the  relationship  between  the  periodicity  and  the  width  of  domain  walls,  particularly  in  PZT 
system.  The  theoretical  snidy  will  also  address  the  grain  size  dependence  of  the  domain  period, 
which  can  shed  some  light  on  the  di^erent  domain  patterns  obtained  from  Transmission  Electron 
Microscopy. 

b)  Study  the  octahedra  tilt  transition  in  the  riiumboheuial  phase  of  PZT  and  the  symmeLy  relation 
of  the  antiphase  boundaries  in  the  R3c  phase  to  interpret  the  TEM  results  reported  in  Appendix  m. 

c)  Exrend  the  model  described  in  Appendix  I  and  n  to  inemporate  the  effect  of  external  electric 
field.  The  change  of  distribution  under  an  electric  field  will  facilitate  the  calculation  on  the  extrinsic 
contributions  from  the  field  induced  phase  transitions. 

d)  Enhance  the  link  between  this  program  and  the  Navy  funded  URI  program  and  try  to  integrate 
the  theoretical  results  with  the  experimental  results.  It  is  expected  that  the  experimental  work  of  the 
URI  program  will  catch  up  in  the  coming  year. 
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A  method  is  introduced  to  determine  the  statistical  divitthitiion  of  enerjjenciitly  Jcsencraie  bttt  geometrically  iti- 
equivalent  states  in  a  temperature  induced  phase  transition  in  solids.  The  method  has  heen  employed  to  calculate  the 
ratio  of  the  rhombohedrai  and  tetragonal  phases  in  the  Pb(Zr,  Ti, . ,  lO,  solid  sohtii><it  of  the  morphotropic  phase  bitund- 
ary  (MPB)  composition.  Our  results  indicate  that  the  MPB  determined  by  Jalle,  CooL  and  Jalfe  i Piezoelectric  Ceramt  i 
(Academic  Press,  London,  1971)  p.  I3A|  from  structural  measurements  should  he  slniied  to  the  rhoinboliedrai  side, 
which  is  more  consistent  with  the  MPB  deterrniried  l•l»n  diclc,.tiic  ,iicasu.ciiiwiiis. 

KEYWORDS:  farroalactric  ceramic.  PZT.  lead  zirconate  tttanate,  iTiorphotropic  phaso  boundary 


A  phase  diagram  representing  the  subsolidus  phase 
relations  in  the  system  PbZrOi-PbTiOj  taken  from  the 
book  by  Jaffe,  Cook  and  JafTe"  is  given  in  Fig.  I.  The 
system  contains  the  well  known  solid  solution  series 
Pb<Zrj,Ti|.,)0]  (PZT)  which  exhibits  a  variety  of  phase 
changes.  In  this  system  the  most  interesting  compositions 
are  for  those  on  or  near  the  so-called  morphotropic 
phase  boundary  (MPB)  (see  the  nearly  vertical  line  on 
Fig.  1).  PZTs  in  this  compositional  region  have  superior 
piezoelectric  properties  and  are  the  primary  materials  cur¬ 
rently  used  in  most  piezoelectric  transducers  and  ac¬ 
tuators.  As  shown  in  Fig.  1  the  ferroelectric  phase  is 
rhombohedrai  on  the  left-hand  side  and  tetragonal  on 
the  right-hand  side  of  the  MPB  line,  respectively  (note: 
there  should  be  a  co-existence  region  of  the  two  phases 
which  is  not  shown  in  this  phase  diagram).  The  MPB  in 
Fig.  1  as  quoted  from  the  book  (JafTe,  Cook  and  JafTe,  p. 
136)  “is  considered  as  that  composition  where  the  two 
phases  are  present  in  equal  quantity”. 

The  MPB  drawn  on  a  phase  diagram  is  a  composi- 


Flg.  I.  PhaM  diagram  of  Pb<Zr,Ti,.,0|)  (after  Jalfe,  Cook  and 
Jeffe,  ref.  II. 


tional  boundary  which  is  defined  as  the  composition  for 
which  the  free  energies  of  two  adjacent  phases  are  equal. 
From  thermodynamics,  the  free  energy  of  a  solid  solu¬ 
tion  system  depends  on  the  following  variables:  lempera- 
ture,  pressure  and  composition.  The  compositional 
variable  in  the  PZT  case  has  a  special  properly,  it  can 
only  be  directly  accessed  in  certain  temperature  range. 
Out  of  that  temperature  range,  the  compositional 
variable  can  not  be  changed  while  keeping  other  variables 
fixed,  as  we  know  that  a  PZT  solid  solution  can  not  be 
formed  at  room  temperature.  In  order  to  change  the  com¬ 
position  at  room  temperature,  one  has  to  first  heat  up  the 
system  so  that  a  single  phase  solid  solution  can  be 
formed,  then  cool  the  system  back  down  to  room  temper¬ 
ature  after  the  new  composition  is  formed.  This  special 
property  of  the  compositional  variable  prevents  direct 
thermodynamic  analysis  with  this  variable,  and  hence  in¬ 
validates  the  transition  hysteresis  argument-'  in  the  ex¬ 
planation  of  the  coexistence  of  two  phases.  Although  the 
phase  transition  fre  n  the  rhombohedrai  structure  to  the 
tetragonal  structure  must  be  of  first  order,  it  can  not  be 
realized  while  keeping  the  temperature  unchanged, 
because  the  composition  is  not  a  directly  accessible  ther¬ 
modynamic  variable  below  a  certain  temperature  limit. 
This  is  to  say  that  we  should  address  the  temperature  in¬ 
duced  transition  process  in  order  to  understand  the 
effects  of  changing  composition  below  u  certain  tempera¬ 
ture  limit. 

There  are  two  types  of  phase  mixing:  one  is  the  mixing 
of  phases  of  dillerent  chemical  compositions  and  the 
other  is  the  mixing  of  phases  with  different  structures  but 
identical  chemical  compositon.  At  a  first  glance,  it  ap- 
(Tcars  that  the  two  problems  seem  to  be  similar,  but  they 
actually  have  quite  different  nature.  The  former  reflects 
the  law  of  mass  conservation  (obeying  the  lever  rule),  but 
the  latter  is  actually  a  statistical  distribution  problem.  It 
is  the  intention  of  this  pajser  to  provide  a  simple  method 
dealing  with  the  latter  case. 

From  thermodynamic  analyses,  it  is  shown  that  tlie 
paraelectric-ferroelectric  transition  for  compositions 
near  the  MPB  is  of  second  order  within  the  experimental 
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error.”  we  therefore  expect  strong  thermal  fluctuations  to 
occur  near  the  transition  temperature  r..  This  thermal 
fluctuation  is  the  sole  driving  force  for  the  phase  transi¬ 
tion.  For  convenience  we  will  work  in  the  order 
parameter  space.  When  the  system  is  cooled  down  near 
Tt  from  the  paraelectric  phase,  the  magnitude  of  the  ther¬ 
mal  induced  instant  polarization  iP,„ul  increases  as  the 
potential  well  around  P=0  becomes  flatter  and  flatter. 
Below  Tc,  the  fluctuation  is  gradually  frozen  and  the 
system  acquires  a  new  configurational  state  with  a  finite 
polarization.  In  a  second  order  phase  transition  the  tran¬ 
sition  process  happens  in  a  continuous  fashion;  but  in  a 
first  order  phase  transition,  a  finite  spontaneous  polariza¬ 
tion  is  obtained  at  the  phase  transition.  The  two  cases  are 
illustrated  in  Fig.  2  for  a  one-dimensional  system,  where 
the  shaded  area  represents  the  level  of  thermal  energy.  In 
a  one-dimensional  system  the  number  of  allowed  (two) 
orientations  of  the  instant  polarization  is  the  same  as  that 
of  the  spontaneous  polarization.  But  in  two-  and  three- 
dimensional  systems,  there  are  no  orientational  con¬ 
straints  for  thermal  fluctuation  (determined  by  the  dimen¬ 
sionality  of  the  problem  only)  although  the  magnitude  of 
Pm  is  regulated  by  the  potential  well  around  the  origin. 
However,  the  allowed  orientations  for  the  spontaneous 
polarization  Pa  are  limited.  The  thick  line  and  the  solid 
dots  in  Fig.  3  illustrated  this  situation  for  a  two-dimen¬ 
sional  problem  representing  a  square  to  rectangular  fer¬ 
roelectric  phase  transition.  There  are  no  orientational 

Tsli4htly>Tc  T  tliehtly  <  Te  T»«llb*lewTc 

First  erder  \  /  \  /  \  / 

\y^y 

Fig.  2.  lUuHraiion  of  the  potential  wdb  and  thermal  fluctuations 
near  the  critical  temperature  r,  for  both  second  and  first  order  phase 
transitions  in  one  dimension.  The  shaded  areas  represent  thermal 
energy. 
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Fig.  3.  Representation  of  fluctuating  state  and  the  ferroelectric  states 
in  order  parameter  space  for  a  two-dimensional  problem,  fa)  The 
thick  and  the  thin  lines  are  the  profiles  for  the  polarization  fluctua¬ 
tions  of  the  square-rectangular  and  square-oblique  phase  transitions, 
respectively,  near  T,.  (b)  The  solid  dots  are  the  degenerated  rec¬ 
tangular  fenodcctric  states  and  the  open  circles  are  the  degenerated 
oblique  ferroeieciric  states,  retpectively,  well  below  T^. 


constraints  for  P.„.<  in  the  fluctuating  state  [profile  is 
represented  by  the  thick  line  in  Fig.  3(a)l  but  only  four 
orientations  for  P.,  (solid  dots  in  Fig.  3(b))  are  allowed 
well  below  T,. 

The  question  vve  try  to  address  in  this  paper  is:  what  is 
the  probability  fi.  for  the  system  in  Fig.  3(a)  to  become 
the  /th  low  symmetry  ferroelectric  states  in  Fig.  3(b)?  For 
the  two-dimensional  problem  mentioned  above,  the 
answer  ^,=  1/4  may  be  obtained  intuitively  since  the  four 
low-temperature  rectangle  states  are  completely 
equivalent,  i.e.,  these  states  are  energetically  degenerate 
and  structurally  identical.  But  in  more  general  situations, 
intuition  fails  to  provide  an  answer,  for  instance,  the  case 
shown  in  Fig.  3  where  the  energetically  degenerated  rec¬ 
tangle  (solid  dots)  and  the  oblique  (open  circle)  states  arc 
structurally  different.  In  this  case  we  need  to  follow  a 
well-defined  mathematical  method  to  calculate  these  pro¬ 
babilities.  Such  a  calculation  is  particularly  useful  for  the 
study  of  PZT  system,  because  at  the  MPB  wc  have  pre¬ 
cisely  such  a  situation,  i.e.,  the  tetragonal  and  rhom- 
bohedral  phases  are  energetically  degenerate. 

it  is  obvious  that  the  probability  calculation  is  only 
meaningful  in  the  fluctuating  state.  After  the  system  be¬ 
ing  frozen  into  one  of  the  low  temperature  states,  ther¬ 
mal  energy  is  no  more  effective  to  carry  the  system  from 
one  configuration  to  another.  Our  calculation  is  based  on 
the  following  assumption; 

The  instant  polarization  P,,„  is  orietuationally 
ergot ic  in  the  fluctuating  state  near  TV. 

In  other  words,  we  assume  that  the  thermal  motion  has 
no  orientational  preference,  although  as  shown  in  Fig.  3, 
the  magnitude  IP,„hI  is  regulated  by  the  potential  con¬ 
figuration  and  is  a  function  of  direction  and  temperature. 
An  immediate  inference  from  the  assumption  is  that  the 
average  total  polarization  of  the  system  is  zero  in  the  fluc¬ 
tuating  state,  or  more  concisely.  <Pi,«>=0  but  <PrnH> 
^0. 

Now  we  try  to  correlate  this  thermal  fluctuation  with 
the  transition  probabilities  to  different  structural  phases. 
The  assumption  tells  us  that  the  trigger  from  the  sur¬ 
rounding  thermal  bath  is  isotropic,  but  the  actual  struc¬ 
tural  change  resulting  from  the  trigger  depends  on  the 
potential  energy  configuration.  Because  only  a  few 
isolated  orientations  for  P,  are  allowed  below  in  the 
order  parameter  space,  we  expect  that  ail  P,n,  oriented  in 
the  vicinity  of  an  allowed  polarization  direction  can 
potentially  develop  into  that  final  polarization  state. 
Hence  we  can  assign  each  allowed  polarization  state  an 
effective  solid  angle  in  a  three-dimensional  order 
parameter  space,  the  probability  pt  of  that  state  being 
formed  under  the  trigger  of  thermal  fluctuation  is 
represented  by  n,/(47r),  where  4it  is  the  normalization 
constant. 

The  next  task  is  to  define  the  boundaries  of  these  solid 
angles,  imagine  we  apply  a  small  dc  field  £  to  lift  the 
degeneracy  of  the  system.  Under  this  field  the  system  wilt 
be  forced  to  develop  into  one  of  the  allowed  low  tempera¬ 
ture  states  whose  polarization  vector  has  the  smallest 
possible  angle  with  the  applied  electric  field.  If  this  field  is 
applied  to  the  system  from  another  angle,  it  can  either  in- 
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duce  the  same  polarization  state  or  a  different  state  de¬ 
pending  on  whether  or  not  the  held  is  still  oriented  inside 
"the  effective  solid  angle  of  that  state.  Considering  at  a 
temperature  slightly  below  7",.  there  are  two  adjacent 
energetically  degenerate  states  represented  by  P,  and  P„ 
respectively,  in  the  order  parameter  space  (P,  and  Pb 
form  an  angle  ^),  the  boundary  of  the  two  effective  solid 
angles  for  these  polarization  states  may  be  defined  by 

£  P,=£  Pb  (1) 

Taking  the  coordinate  systems  as  shown  in  Fig.  4  we  can 
write  eq.  (I)  in  the  following  form. 

P.COS^*PbCOS(<fr--^).  (2) 


Interestingly,  the  magnitude  of  the  electric  field  has  drop¬ 
ped  out  in  eq.  (2),  hence  if  we  take  the  limit  £->0,  the 
result  in  eq.  (2)  will  still  hold. 

In  reality,  the  problems  of  interest  are  those  cases  for 
which  P,“Pb,  therefore,  eq.  (2)  can  be  further  simplified 
to 

*»«/2.  (3) 


Equation  (3)  states  that  the  boundary  of  two  effective 
solid  angles  is  a  plane  in  the  order  parameter  space  pass¬ 
ing  through  the  bisector  of  the  angle  ^  and  perpendicular 
to  the  plane  containing  P.  and  Pk.  If  there  exist  more 
than  two  degenerate  states,  a  boundary  may  be  defined 
for  each  pair  of  adjacent  effective  solid  angles  using  eqs. 
(2)  or  (3). 

After  the  boundaries  are  defined,  the  remaining  task 
becomes  straightforward.  We  draw  a  polyhedron  in  the 
order  parameter  space  surrounding  the  origin,  whose 
edges  are  on  the  solid  angle  boundaries  defined  by  eq.  (3). 
For  a  PZT  of  MPB  composition,  the  degenerated  fer¬ 
roelectric  states  include  both  rhombohedral  and 
tetragonal  phases.  Assuming  for  the  simplest  case  that 
Pjmp^mg,  then  the  corresponding  polyh^ron  is  shown 
in  Fig.  5.  The  solid  angle  calculation  can  be  written  in 
terms  of  a  surface  integral  on  each  of  the  corresponding 
surface  which  subtends  that  solid  angle,  for  example 
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rig.  4.  Coordtaautytttmiistdi«tiitcalculaiioiuorcqs-(2)aiid(3). 


Fi|.  5.  The  crobabiliiy  diMribuiian  potyhedron  for  PZT  of  (he  MPB 
compotiiion,  OT ^OR.  The  solid  engic  subtended  by  the  surface 
IIKH  with  respect  to  the  origin  O  divided  by  4]i  represents  (he 
probability  of  the  system  being  transformed  into  (00/*„)  state. 


where  dr  is  the  area  element  on  the  surface  OHIJK  and  r 
is  the  distance  between  dr  and  the  center  point  T  on  that 
surface  (Fig.  S).  The  integration  of  eq.  (S)  can  be  carried 
out  to  give  an  analytic  solution 


00111(11=^4  arcsin 


(6) 


From  this  result  and  the  arguments  given  above  we  con¬ 
clude  that  for  PZT  of  the  MPB  composition,  the 
probability  ratio  for  a  system  to  be  frozen  into  rhom- 
bohcdral  and  tetragonal  phases  from  the  fluctuating  state 
is  given  by 

.  /3- V 

Rhombohedral  1 -6QuHiJK/4)t  *  *fcstn  ^  ^ 

■  ,  II  ■  _  —.—  as  i-  I  .■■a  — gi- 

Tetragonal  6Qoiiijk/4*  _  .  /3-*3 

fiarcsin  I — g — 
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A  ceramic  system  may  be  treated  as  an  ensemble  of  do¬ 
mains,  and  each  domain  can  be  considered  to  be  the 
system  we  have  discussed  above.  Then,  eq.  (7)  represents 
the  molar  ratio  of  the  rhombohedral  and  tetragonal 
phases  for  a  PZT  ceramic  of  the  MPB  composition  l  his 
molar  ratio  can  be  calculated  directly  from  X-ray  ditlrac- 
tion  intensities,  therefore  eq.  (7)  provides  a  criterion  for 
the  determination  of  the  MPB. 

Thermodynamics  tells  us  that  maximum  values  of 
many  physical  quantities  should  appear  at  the  MPB  com¬ 
position  due  to  the  existence  of  maximum  number  of 
energetically  degenerated  states.  However,  several  ex¬ 
perimental  results  reveal  that  these  maxima  do  not  match 
with  the  MPB  on  Fig.  I  but  often  slightly  deviate  to  the 
rhombohedral  side.^'^’  Our  analysis  above  gives  a 
reasonable  explanation  to  this  controversy.  Naturally,  if 
we  use  the  ratio  of  3:2  instead  of  1:1  as  the  criterion  for 
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the  MPB,  the  MPB  line  on  Fig.  I  would  shift  to  the  rhom- 
bohedrai  side. 

The  novel  idea  presented  in  this  paper  is  to  take  into  ac¬ 
count  the  geometrical  constraints  in  the  calculation  of 
the  statistical  distributions  of  those  energetically 
degenerated  states  in  solids.  Through  this  paper  we  also 
intend  to  bring  people’s  attention  on  treating  solid 
systems,  those  thermodynamic  theories  developed  for  gas 
and  liquid  systems  may  not  be  valid  and  should  be 
augmented  to  incorporate  the  characteristics  of  solids. 

Finally,  we  like  to  point  out  that  a  special  property  per¬ 
taining  to  this  problem  has  been  used  implicitly,  i.e.,  the 
homotopy  mapping  between  the  real  space  and  the  order 
parameter  space.  Strictly  speaking,  thermal  motion  oc¬ 
curs  in  real  space  not  in  order  parameter  space.  But 
because  the  order  parameter  is  a  vector  in  our  problem, 
there  is  a  one  point  to  one  point  mapping  between  the 


real  space  and  the  order  parameter  space,  therefore  we 
did  not  emphasize  the  difference  between  the  two  spaces 
in  the  text.  However,  if  the  order  parameter  is  not  a  vec¬ 
tor,  one  can  not  use  the  theory  developed  in  this  paper. 
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A  statistical  model  is  proposed  to  address  the  problem  of  two-phase  coexistence  near  the  morphotro¬ 
pic  phase  boundary  (MPB)  in  PbiZri  ..Ti.  )0)  solid-solution  series.  Functional  forms  for  the  molar  frac¬ 
tions  of  tetragonal  and  rhombohedral  phases  inside  the  coexistence  region  are  obtained,  which  may  be 
used  to  replace  the  lever  rule  to  describe  the  phase  mixing  in  a  complete  binary  solid-solution  series 
without  solubility  gap.  The  model  predicts  that  the  width  of  this  coexistence  region  is  inversely  propor¬ 
tional  to  the  volume  of  each  element  in  the  sutistical  ensemble.  In  addition,  the  shift  of  the  MPB  com¬ 
position  from  the  composition  of  equal  molar  fraction  of  the  two  coexisting  phases  is  found  to  be  propor¬ 
tional  to  the  width  of  the  coexistence  region.  Several  existing  controversial  experimental  observations 
can  be  reconciled  by  this  model. 


L  INTROOUenON 

The  most  widely  used  piezoelectric  certmic  todsy  is 
lesd  zirconate  titanate  (PCT),  a  solid  solution  of  PbZtOj- 
PbTiO},  with  compositions  near  the  morphotropic  phax 
boundary  (MPB).  A  MPB  is  defined  as  a  compositional 
phase  bMndary  at  which  the  two  adjacent  phases  in  a 
phase  diagram  have  equal  Gibbs  free  energy.  The  phase 
diagram  determined  by  Jaflb,  Cook,  and  Jaflb'  from  x- 
ray-difl^action  measurements  is  shown  in  Fig.  1.  The 
MPB  on  this  phase  diagram  was  considered  to  be  the 
composition  at  which  the  amount  of  tetragonal  and 
rhombohedral  phases  is  equal.'  However,  it  was  pointed 
out^  that  the  molar  ratio  of  the  two  low-temperature 
phases,  i.e.,  rhombohedralitetngtuul,  should  be  3:2  in¬ 
stead  of  1:1  at  the  MPB  (defined  by  equating  the  free  en¬ 
ergies  of  the  two  leases),  which  provides  an  explaiution 
for  the  discrepancy  between  the  MPBs  determined  by  us¬ 
ing  dielectric  maximum  and  from  x-ray-difliraction  inten¬ 
sities. 

Historically,  the  exact  composition  of  MPB  in  PZT  ha* 
never  been  precisely  defined;  it  ranges  from  4S-S0  mol  % 
of  PbTiOy.^'*  Ttere  is  a  coexistence  region  of  the 
tetragonal  and  rhombohedral  phases  whose  width  is  also 
not  well  defined,*'*  ranging  from  2-  IS  mol  %  of  PbTiOj 
around  the  composition  Ti/Zr*4g/32. 

From  many  years  of  study  on  the  PZT  system,  the  fol¬ 
lowing  two  factt  are  well  accepted: 

(1)  The  PZT  system  is  a  complete  binary  solid  solution 
of  PbTtOy  and  PbZt03  without  solubility  gaps. 


(2)  Below  the  paraekctric-ferroelectric  transition  tem¬ 
perature  there  exisu  a  coexistence  region  of  the  tetrago¬ 
nal  and  rhombohedral  phases  near  the  MPB  composition, 
although  the  width  of  this  region  is  still  a  debatable  issue. 
Adding  small  amounts  of  dopant  can  shift  the  .VfPB  and 
increase  the  width  of  the  coexistence  region. 

The  lever  rule,  obtained  from  mass  conservation,  has 
been  used  to  describe  fact  (2)  above.^’*  Although  the  data 


FIO.  1.  Phase  diagram  of  PbZiOj-PbTiO,  P/I  '•  I'J- 
solution  series  by  Jallie,  Cook,  and  Jaffe  (Ref.  1). 
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fitting  appears  to  be  reasonable,  the  two  edge  composi- 
rtons  X ,  and  in  the  lever  rule  actually  specify  a  solubil* 
ity  gap,  which  is  in  contradiction  with  fact  (1)  above.  An 
alternative  explanation  of  the  phase  coexistence  [fact  (2) 
aSove]  is  to  use  the  transition  hysteresis  argument  within 
the  context  of  a  first-order  phase  transition  by  taking  the 
composition  as  an  independent  variable. This  is  also 
inadequate  because  the  composition  variable  is  already 
frozen  in  the  temperature  region  (T  <375*0  where  the 
two  ferroelectric  phases  exist;  the  system  only  experi¬ 
ences  a  diffusionless  structural  phase  transition  at  the  Cu¬ 
rie  point.  Solid-state  reaction,  which  is  needed  to  change 
the  composition  of  PZT,  cannot  occur  until  above  800  *C. 
In  other  words,  once  a  solid  solution  is  formed  at  high 
temperature  (r>  800*0,  the  chemical  composition  can¬ 
not  be  changed  at  low  temperatures,  but  the  system  can 
have  temperature-induced  diffusionless  structural  phase 
transitions.  Because  the  compositional  degree  of  freiKlom 
has  been  frozen  in  the  temperature  region  of  interest,  the 
transition  hysteresis  concept  is  not  feasible.  In  addition, 
both  explanations  mentioned  above  lead  to  a  definite 
width  of  the  coexistence  region  which  has  not  been  ob¬ 
served  experimentally. 

Many  physical  properties  of  PZT  reach  their  max¬ 
imum  or  minimum  at  the  MPB,  which  is  also  not  weU  un¬ 
derstood.  From  the  definition  of  MPB,  the  two  low- 
temperature  structural  phases  are  energetically  degen¬ 
erate  at  the  MPB  composition,  and  it  is  conceivable  that 
the  electric  field  or  stress-driven  phase  transitions  be¬ 
tween  the  two  ferroelectric  structures  are  possible  for  the 
PZT  of  composition  near  the  MPB.  This  field-induced 
phase  transition  could  contribute  substantially  to  the  ob¬ 
served  phenomena.  In  order  to  quantify  this  contribu¬ 
tion.  one  must  know  the  exact  fraction  of  the  two- 
structural  phases  for  a  given  composition  in  the  low- 
temperature  regime.  Looking  at  Clascal  thermodynam¬ 
ic  and  sutistical  theories,  one  finds  that  they  cannot  be 
directly  applied  to  address  our  problem  because  the  MPB 
is  defined  by  equating  the  fne  energies  of  the  tetragonal 
and  rhombohedral  phases,  hence  the  energy  difference 
(which  is  the  only  criterion  in  classical  statistics)  is  zero. 
On  the  other  hand,  after  a  low-temperature  structure  is 
formed  from  the  paraelectric-ferroelectric  phase  transi¬ 
tion,  the  system  may  be  "locked**  into  this  structure  be¬ 
cause  of  the  existence  of  a  transitioo  barrier  between  the 
two  low-temperature  phases  (this  is  reflected  in  Fig.  1  as 
the  nearly  vertical  line  of  the  MPB).  This  “locked  struc¬ 
ture”  is  thermodynamically  metastable  below  the  transi¬ 
tion  temperature.  Therefore,  in  reality  we  are  not  exactly 
dealing  with  the  absolute  thermodynamic  equilibrium  at 
temperatures  well  below  Tf.  It  is  our  opinion  that  the 
coexistence  of  the  two  low-temperature  phases  in  the 
PZT  system  is  a  result  of  quenched-in  thermal  fluctua¬ 
tions.  Following  this  idea,  the  two  observations  in  the 
PZT  system  mentioned  above  can  be  explained  satisfacto¬ 
rily,  a^  the  fractions  of  the  two  coexisting  phases  can  be 
quantifled. 

n.  THE  MODEL 

As  a  starting  point,  we  assume  that  the  partitioning  of 
the  two  phases  does  not  change  with  temperature  after 


being  determined  at  the  paraelectnc-ferroelectnc  phase 
transition.  We  therefore  only  need  to  calculate  the  parti¬ 
tioning  near  T^- 

The  question  we  are  trying  to  address  here  is  really  the 
accessibility  of  all  the  low-temperature  states  during  the 
phase  transition.  If  we  assume  that  all  14  states  (8  in  the 
rhombohedral  phase  and  6  in  the  tetragonal  phase)  are 
identical,  then  the  molar  ratio,  rhombohedral:tetragonal. 
should  be  4:3  at  the  MPB.  But  obviously  these  14  states 
are  not  identical,  the  two  low-temperature  structures  are 
geometrically  inequivalent.  In  Ref.  2  we  have  applied 
this  geometrical  constraint  in  the  statistical  calculation  of 
this  ratio  at  the  MPB,  which  is  close  to  3:2,  and  intro¬ 
duced  the  concept  of  probability  polyhedron  for  systems 
with  a  vector  order  parameter.  This  concept  may  be  used 
to  calculate  the  distribution  of  energetically  aegenerate 
but  geometrically  inequivalent  states  resulting  in  a 
second-order  phase  transition  (here  the  transition  refers 
to  the  paraelectric-ferroelectric  transition  but  not  the  one 
between  the  Iwo  low-temperature  phases).  In  this  paper 
we  try  to  extend  the  model  discuss^  in  Ref.  2  to  account 
for  the  coexistence  of  two  energetically  nondegenerate 
and  geometrically  inequivalent  phases  by  incorporating 
the  classical  statistics. 

In  order  to  understand  the  underlying  physics  of  the 
present  problem,  we  first  study  the  driving  force  for  the 
phase  transition.  It  has  been  shown  using  phenomeno¬ 
logical  theory  that  the  paraelectric  to  ferroelectric  phase 
transition  in  PZT  for  compositions  in  the  vicinity  of  the 
MPB  is  of  second  order."  Therefore,  thermal  fluctua¬ 
tions  are  the  driving  force  for  this  phase  transition.  Inev¬ 
itably,  these  fluctuations  will  also  play  a  key  role  in  the 
prolmbility  distribution  of  the  tetragonal  and  rhom¬ 
bohedral  phases  during  cooling  through  the  Curie  point. 

A.  The  prebablUty  polyhedroa 

The  construction  of  probability  polyhedron  was  de¬ 
scribed  in  Ref.  2.  The  only  assumption  made  there  was 
the  orientational  ergodicity  of  thermal  fluctuations  which 
is  valid  for  either  long-time  observation  of  a  single  system 
or  instant  observation  of  a  statistical  ensemble. 

In  order  to  generalize  the  idea  of  probability  polyhed¬ 
ron  we  give  an  equivalent  definition  below.  Considering 
the  fact  that  the  number  of  surface  planes  in  the  polyhed¬ 
ron  is  equal  to  the  number  of  allowed  polarization  direc¬ 
tions,  it  is  equivalent  to  say  that  in  the  fluctuating  sUte 
each  of  the  polarization  sutes  occupies  an  effective  solid 
angle  in  the  order  parameter  space,  which  is  equal  to  the 
solid  angle  subtended  by  the  polyhedron  surface  who* 
plane  normal  coincides  with  that  polarization  directive. 
The  probability  of  attaining  this  polarization  state 
cooling  throu^  the  phase  transition  is  equal  to  its 
effective  solid  angle  divided  by  4ir,  the  normalization  con¬ 
stant. 

Now  we  use  this  concept  to  describe  the  distribution  of 
polarization  stttes  for  the  PZT  of  compositions 
sufficiently  far  from  the  MPB.  In  this  case,  the  low- 
temperature  phase  is  either  tetragonal  or  rhombohedral. 
On  the  tetragonal  side  (Ti  rich)  of  the  MPB,  the  probabil¬ 
ity  polyhedron  is  a  cube  with  each  of  the  six  variants. 
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(ip.0,0).  (O.ip.O),  and  (0,0, ±p),  occupying  an  efTective 
solid  angle  of  O,  =  2ir/3  (i  =  1  -6),  where  p  is  the  magni¬ 
tude  of  the  polarization.  Hence,  the  fraction  being 
transformed  into  each  of  the  polarization  states  in  an  en¬ 
semble  (a  ceramic  can  be  treated  as  an  ensemble  of  grains) 
is  equal  to  fl,  /( 4ir )  =  I  /6,  and  the  total  polarization  of 
the  ensemble  is  zero  in  the  order  parameter  space.  Be¬ 
cause  of  the  one  to  one  point  mapping  between  the  order 
parameter  space  and  the  real  space,  the  total  polarization 
in  the  real  space  is  also  zero  for  a  statistical  ensemble  of 
particles.  The  same  can  also  be  said  for  the  compositions 
on  the  rhombohedral  side  (Zr  rich)  of  the  MPB,  for 
which  the  probability  polyhedron  is  an  octahedron.  The 
effective  solid  angle  for  each  variant  is  tr/2  and  the  prob¬ 
ability  for  each  polarization  state  is  (Tr/2)/(4rr)=‘l/8. 
Again,  the  net  polarization  ot  the  ensemble  of  rhom¬ 
bohedral  phase  particles  is  zero.  When  the  PZT  compo¬ 
sition  is  close  to  the  MPB  composition,  the  two  low- 
temperature  structural  phases  can  coexist  and  the  proba¬ 
bility  polyhedron  will  have  14  faces.  In  general,  the 
tetragonal  states  and  the  rhombohedral  states  are  not  de¬ 
generate  unless  the  composition  is  exactly  on  the  MPB; 
therefore,  we  expect  the  effective  solid  angle  representing 
the  probability  of  each  polarization  state  to  change  with 
composition. 

There  is  an  energy  difference  between  the  two  structur¬ 
al  phases  when  the  composition  is  not  exactly  on  the 
MPB;  we  introduce  an  anisotropy  factor  6  to  describe 
this  situation.  As  shown  in  Fig.  2,  the  distance  of  the  sur¬ 
faces  (corresponding  to  different  phases)  from  the  center 
point  of  the  polyhedron  is  represented  by  r,-  (=07  or 
OR),  which  controb  the  solid  angle  subtended  by  the  sur¬ 
face.  In  other  words,  determines  the  probability  of  at¬ 
taining  a  specific  polarization  sute  using  the  concept  we 
have  introduced  above;  r,  must  be  a  function  of  energy. 
At  the  MPB  composition,  the  tetragonal  and  the  rhom¬ 
bohedral  phases  are  degenerate,  and  we  have  rj^r^.  In 
thu  case,  the  actual  magnitude  of  r,  does  not  matter  since 
the  solid  angle  partitioning  u  independent  of  r,  .  But  for 
the  nondegenerate  case,  rj  and  are  different  and 
should  depend  on  the  energies  of  the  two  low- 


FIO.  2.  Probability  polyhedron  for  PZT  systein  with  the  pa- 
rameter  8  in  the  range  1  — l/v'3>8>  I— 2/v^3.  At  MPB, 
Or=OR  and  6=0. 


temperature  phases.  We  define  the  distnbution  anisotro¬ 
py  parameter  6  as  follows: 


'■r 


which  is  a  function  of  composition  only  when  tempera¬ 
ture  and  pressure  are  fixed. 

S=0  represents  the  MPB  composition,  at  which  the 
probability  ratio  of  the  rhombohedral  and  tetragonal 
phases  is* 

fn  _  ff— 6arcsin[(3- V' 3)/6l 
fr  6arcsin[(3  — v'li/b] 

where  /j  and  fj  are  the  probabilities  of  the  rhom¬ 
bohedral  and  tetragonal  phases.  This  ratio  is  close  to  3:2 
instead  of  1:1  given  by  Jaffe,  Cook,  and  Jaffe.' 


B.  The  distribntfon  functions  for  the  coexisting  phases 
in  PZT  of  compoaition  near  the  MPB 

There  are  upper  and  lower  bounds  for  the  value  of  6. 
When  8  decreases,  the  probability  of  transforming  in^o 
the  rhombohedral  phase  also  decreases  until  5  =  1  -  v  3; 
for  S  less  than  thb  critical  value,  the  low-temperature 
phase  can  only  be  tetragonal  because  the  polyhedron  be¬ 
comes  a  cube.  On  the  other  hand,  when  5  increases,  the 
probability  of  transforming  into  the  tetragonal  phase  de¬ 
creases,  the  upper  limit  for  6  is  I  — l/v'l.  for 
6^1  —  1  /v^;  the  polyhedron  becomes  an  octahedron, 
which  means  that  the  system  can  only  be  rhomboh^ral. 
There  is  another  special  value  of  6,  6=l-2/v^3.  at 
which  the  representative  surface  of  the  rhombohedral 
phase  on  the  probability  polyhedron  changes  its  shape 
from  a  six-sitM  polygon  to  a  right  triangle,  and  the  j 
representative  surfiice  for  the  tetragonal  phase  changes 
from  a  square  to  an  eight-sided  polygon.  Therefore,  in 
calculating  the  effective  solid  angte  for  each  of  the  polar¬ 
ization  states,  one  must  use  Fip.  2  and  3  for  the  cases  of 
6>1— 2/Vl  and  6<1— 2/V3,  respectively.  It  can  be 


FIO.  3.  Probability  polyhedron  for  PZTjystem  wnti  ihe  pa¬ 
rameter  6  in  the  range  1— 2/v'3>8>  l  —  y'Xrr^r,,. 
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C.  Tht  physical  mcaaiag  of  6  aad  the  width 
of  the  cocsiateacc  regioa 


In  the  PZT  system,  the  free  energies  of  the  tetragonal 
and  rhombohednl  phases,  Gj-  and  G^,  depend  monotoni- 
cally  on  composition.  The  two  free  energies  cross  each 
other  at  the  MPB.'^  From  thermodynamics,  at  equilibri¬ 
um  only  one  of  the  two  low-temperature  phases  is  stable 
except  at  the  MPB  composition.  However,  in  the  vicinity 
of  the  MPB,  one  expects  the  energy  difference  of  the  two 
low-temperature  phases  to  be  very  small;  thermal  energy 
can  intr^uce  some  uncertainties  in  the  distribution  /,, 
which  obeys  the  canonical  distribution.  Since  we  have  as¬ 
signed  the  solid  angle  n,  to  be  proportional  to  the  distri¬ 
bution  function  /,  ,  11^  nuy  be  written  as  follows: 


n<«/.ccexp 


kT, 


(3) 


where  le  is  Boltzmann’s  constant,  G<  and  Gg  are  the  sys¬ 
tem  free  energies  for  the  /tb  low-temperature  structural 
phase  and  for  the  lowest  energy  phase  (“ground  state”), 
respectively.  For  a  given  surface  area,  the  solid  angle  it 
suMends  with  respect  to  a  given  point  in  space  is  inverse¬ 
ly  proportional  to  the  square  of  the  distance  between  the 
surface  and  that  point,  i.e.,  12, « 1  /rf.  Hence,  firom  Eq. 
(S)  the  distance  variable  r,  can  be  written  in  terms  of  the 
free-energy  difference*. 


r<« 


3/n: 


-exp 


2*r, 


(6) 


By  substituting  Eq.  (6)  into  Eq.  (1)  the  parameter  6  be¬ 
comes 


5-1— exp 


Gk-Gt 

IkTg 


(7) 


Using  Eq.  (7)  and  the  limiting  values  of  5,  one  can  calcu¬ 
late  the  required  energy  difference  in  order  to  obtain  a 
single  phase  sute: 

|G*~G7-|>ifer,,ln3 .  (8) 


Therefore,  the  width  of  the  coexistence  region  depends  on 
how  fast  the  energy  dUference  Gr~Gg  changes  with 
composition.  An  important  conclusion  can  be  drawn 
from  Eq.  (8);  the  width  of  the  MPB  depends  on  the 
votunm  of  the  element  in  the  ensemble  (for  example,  the 
particle  size  in  a  powder  system).  This  is  because  the 
ftee-energy  diflhrenM  on  the  left-hand  side  of  Eq.  (8)  is  an 


I 

extensive  variable  while  the  product  on  the  right-hand 
side  is  an  intensive  variable,  lliis  point  can  become  more 
transparent  if  we  recast  Eq.  (8)  in  the  following  form: 

=  ,  (9) 

where  gr  and  gg  are  the  free-energy  densities  of  the  two 
phases  and  v  denotes  the  volume  of  each  element  in  the 
statistical  ensemble.  Because  the  free-energy  difference  is 
small  for  PZT  compositions  sufficiently  close  to  the  MPB, 
we  may  write  the  free-energy  density  difference  in  terms 
of  a  series  expansion  around  the  MPB  composition: 

<1*1 

where  x  is  the  composition  variable  and  Xq  is  the  MPB 
composition,  and 

Note  that  Eq.  (10)  is  a  mathematical  representttion  but 
not  the  Landau  free  energy,  and  there  are  no  symmetry 
'constraints  for  the  expansion  coefficients. 

In  a  linear  approximatioa,  i.e.,  talcing  a,  —0  for  a#1, 
the  width  of  MPB  5x  can  be  derived  using  Eqs.  (9)  and 
(10): 


Ax- - ^ln3.  (12) 

a,v 

Equation  (12)  indicates  to  us  that  the  width  of  the  coex¬ 
istence  region.  Ax,  is  inversely  proportional  to  the 
volume  of  the  sutistical  element.  Suppose  Ax  is  0.1  for  a 
particle  size  of  0.1  ftm,  then  it  would  be  0.0123  for  a  par¬ 
ticle  size  of  0.2  itm.  Therefore,  it  is  not  surpri^  that 
the  values  of  Ax  obtained  by  diffinent  processiog  tech¬ 
niques  are  quite  different.  It  is  also  conceivable  that  for  a 
well-sintered  ceramic  system  or  a  single-crystal  system. 
Ax  will  be  too  shmII  to  be  detectaMe  with  the  available 
experimental  techniques,  which  gives  an  explanation  as 
to  why  the  coexistence  could  not  be  observed  in  some  ex¬ 
periments,  especially  in  a  single-crystal  system. 

III.  COMPARISON  WITH  EXPERIMENTS 

Using  the  linear  approximation,  we  can  rewrite  Eq.  (7) 
in  the  following  form: 


I 
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6*  I  -exp 


(x  — Xr,)ln3 

Ix 


(13) 


Here  Xg  is  the  MPB  composition  and  Ax  is  the  width  of 
the  coexistence  region  as  mentioned  above.  These  two 
parameters  can  be  obtained  from  fitting  experimental 
data  to  the  distribution  functions  given  by  Eqs.  (3)  and 
(4).  In  practice,  the  molar  ratio  of  the  two  coexisting 
phases  is  measured  from  the  integrated  intensities  of  the 
x-ray-diffraction  peaks  of  the  rhombohedral  and  tetrago¬ 
nal  phases.  The  value  of  Xg  may  be  obtained  using  Eq. 
(2),  viz.,  at  X  =Xg  the  intensity  ratio  of  the  rhombohedral 
phase  to  the  tetragonal  phase  is  roughly  3;2.  Because  the 
x-ray  technique  has  a  limited  accuracy  for  a  second  phase 
of  less  than  a  few  percent,  especially  when  the  diffraction 
peaks  are  not  well  separated,  it  is  difficult  to  measure  the 
coexistence  width  Ax,  In  order  to  overcome  this 
difficulty,  a  useful  relation  is  given  below; 
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FIG.  4.  The  molar  fractions  of  the  rhombohedral  and  tetrag¬ 
onal  phases  inside  the  coexistence  compositional  region  for 
0. lPho_^K.o.  il2n| /jhibj/j )02. 9s“’i0.9— X  iPbZrO)— xPbTiOi  solid 
solution.  The  squares  and  circles  are  the  experimental  data 
from  Ref.  8,  the  solid  curves  were  fitted  from  the  lever  rule,  and 
the  dotted  curves  were  calculated  using  the  proposed  model  un¬ 
der  linear  approximation. 


w  0.053 Ax  ,  (14) 

where  x„  is  the  equal  fraction  composition  at  which  the 
molar  ratio  of  the  two  phases  is  1:1,  which  can  be  easily 
obtained  from  x-ray-diffraction  measurements.  Equation 
(14)  was  obtained  by  substituting  Eq.  (13)  into  Eq.  (3)  and 

letting /r“  1/2" 

For  the  pure  PZT  system,  the  coexistence  region  is  not 
only  very  narrow  but  also  very  sensitive  to  the  processing 
procedures:  there  are  no  reliable  diffraction  data  in  the 
literature.  But  when  the  system  is  slightly  doped,  the 
coexistence  region  becomes  wider  and  the  peaks  become 
easily  identifiable.  As  an  example,  we  examine  the  exper¬ 
imental  data  of  Hanh,  Uchino,  and  Nomura,*^  which  is 
for  the  solid-solution  system 

0.  IPboi  9K0,  |(Zni/3Nbj/j)02.9j 

-(0.9— x)PbZi03-xPbTi0, . 

The  squares  and  circles  in  Fig.  4  are  the  measured  molar 
fractions  of  the  rhombohedral  and  the  tetragonal  phases, 
respectively;  the  solid  curves  were  obtained  by  Hanh, 
Uchino,  and  Nomura  firom  fitting  the  experimental  data 
to  the  lever  rule,  and  the  dotted  curves  are  from  the 
current  model.  Linear  approximation  [Eq.  (13)}  was  used 
in  the  calculation  and  the  two  parameters  Xg  =‘0.5027 
and  Ax  =>0.2066  were  fitted  to  the  experimental  data  us¬ 
ing  the  nonlinear  Levenberg-Marguardt  method.  The 
two  special  compositions  Xg  (MPB)  and  x„  are  also  given 
in  Hg.  4  as  references. 

Generally  speaking,  linear  approximation  is  valid  only 
when  the  two  free-energy  curves  are  relatively  straight  as 
a  function  of  composition  near  the  crossover  point  of  the 
two  free  energies.  It  is  expected  that  the  calculated  coex¬ 
istence  region  could  become  slightly  narrower  when  the 
full  expansion  in  Eq.  (10)  is  used. 

From  Fig.  4  one  may  find  that  the  lever  rule  seems  to 
give  a  good  fit  to  the  limited  experimental  data  points, 
however,  the  derivatives  of  the  distributions,  d/ji  /dx  and 


dff/dx  are  discontinuous  at  the  compositions  x  j  and  x., 
which  represent  the  existence  of  a  solubility  gap  between 
x,  and  Xj.  This  is  in  contradiction  with  the  complete 
solubility  of  the  system.  On  the  other  hand,  our  model 
not  only  provides  excellent  fit  to  the  experimental  data, 
but  also  eliminates  such  derivative  discontinuities,  whicli 
makes  it  more  suiuble  for  describing  the  phase  mixing  in 
complete  solid-solution  systems. 

In  order  to  further  illustrate  the  difference  between  the 
two  theories,  let  us  look  at  a  binary  system  AC-BC  and 
assume  they  form  solid-solutions  a  and  B  for  .4 -rich  and 
ff-rich  compounds,  respectively.  Then,  for  any  given 
composition  x  inside  the  coexistence  region  of  a  and  /?, 
we  have  the  two  theories  describe  the  following  situa¬ 
tions: 

(a)  Lever  rule 

xAC+ll-x)BC  =/„ ,4, , C( a  structure ) 

-*2^0  structure)  ; 

(15) 

(b)  Present  model 

xAC+ll-x)BC==/aAj,Bi-,C(a  structure) 

CO  structure) 

(16) 

The  lever  rule  specifies  a  solubility  gap  from  x,  to  x.,  and 
the  two  coexisting  phases  have  different  chemical  compo¬ 
sitions  as  shown  on  the  RHS  of  Eq.  (15).  On  the  other 
hind,  our  model  was  derived  from  the  complete  solubility 
of  AC  and  BC,  indicated  on  the  RHS  of  Eq.  (16).  and  the 
phase  coexistence  was  considered  as  a  frozen-in  second 
metastable  phase  from  thermal  fluctuations. 

IV.  SUMMARY  AND  CONCLUSIONS 

A  theoretical  treatment  is  proposed  to  calculate  t lie 
molar  fractions  of  the  rhombohedral  and  tetragonal 
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phases  near  the  MPB  in  a  PZT  system.  Under  the  as- 
Aimption  that  the  partitioning  of  the  two  low- 
temperature  phases  is  determined  at  the  paraelectric- 
■*  ferroelectric  transition,  the  distribution  functions  are  re¬ 
lated  to  the  effective  solid  angles  associated  with  the  low- 
temperature  phases  in  the  order  parameter  space.  Ana¬ 
lytic  forms  were  obtained  for  the  molar  fractions  of  the 
two  low-temperature  phases  inside  the  coexistence  re¬ 
gion.  These  molar  fractions  depend  on  a  single  parame¬ 
ter  6,  which  is  a  function  of  the  free-energy  difference  of 
the  two  low-temperature  phases.  Besides  the  energy  con¬ 
siderations,  the  geometrical  constraints  of  a  solid  struc¬ 
ture  have  been  incorporated  in  the  statistical  calculations. 

Using  this  model,  the  two  contradictory  facts  of  the 
PZT  system  mentioned  in  the  introduction  can  be  recon¬ 
ciled,  as  the  cut  rent  model  allows  for  the  phase  coex¬ 
istence  in  a  complete  solid-solution  system  without  solu¬ 
bility  gaps.  In  addition,  the  controversy  about  the 
undefined  width  of  the  coexistence  region  may  also  be  ex¬ 
plained  using  Eq.  (9),  which  states  that  the  width  of  the 
coexistence  region  is  inversely  proportional  to  the  volume 
of  the  element  in  a  statistical  ensemble  (such  as  the  parti¬ 
cle  size  in  a  powder  system).  Because  the  particle  (or 
grain)  size  depends  strongly  on  the  processing  technique. 


a  unique  value  for  the  width  Xx  cannot  be  well  defined. 
The  coexistence  should  not  occur  in  a  single-crystal  sys¬ 
tem,  which  is  consistent  with  the  experimental  observa¬ 
tions. 

The  model  was  applied  to  the  experimental  data  of 
Hanh,  Uchino,  and  Nomura, ‘  which  is  for  the  complete 
solid-solution  system,  0.1Pbo,Ko  |(Zn,/3Nb2/3)02  ,5 
-(0.9— X  )PbZr03-xPbTi03,  and  compared  with  the 
fitting  using  the  lever  rule.  Although  both  theories  pro¬ 
vide  good  fit  to  the  experimental  data,  the  solubility  gap 
specified  by  the  lever  rule  makes  it  unsuitable  to  this 
problem,  while  the  current  proposed  treatment  can  elimi¬ 
nate  the  two  unphysical  kinks  in  the  distribution  func¬ 
tions  at  X,  and  Xj  given  by  the  lever  rule.  Therefore  the 
proposed  model  is  more  consistent  with  the  nature  of 
complete  soiid-solution  systems.  Several  predictions  were 
made  from  the  current  model,  including  the  relationship 
between  the  width  Ax  and  the  particle  volume,  which 
await  further  experimental  verifications. 
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Abstract 

A  transmission  electron  microscopy  (TEM)  study  was  performed  on  Pb(Zi.Ti)03 
compositions  within  the  R3m  R3c  phase  region.  The  low  temperature  phase  is  owing  to  a 
displacive  phase  transition  involving  oxygen  octahedral  tilts.  The  associated  superlattice  reflections 
are  detectable  by  electron  diffraction.  Dark  Field  diffraction  contrast  imaging  of  the  superlattice 
reflections  reveals  antiphase  boundaries  associated  with  octahedral  tilt  domains.  Interaction 
between  the  octahedral  tilt  antiphase  boundaries  and  the  ferroelectric  domain  structures  of  the  R3c 
phase  is  studied  and  discussed . 


L 


The  perovskite  solid-solution  between  end-member  PbTlOj  and  PbZr03  is  the  basis  of 
important  technological  ceramics  used  in  piezoelectric,  pyroelectric  and  electro-optic 
applications.  1*3  The  phase  diagram  of  PbZr03-PbTi03  is  illustrated  in  Figure  1.  The  phase 
diagram  contains  a  variety  of  displacive  |riiase  transitions  and  there  are  antiferroelectric  and  various 
ferroelectric  phases  in  the  low  teiiq)erBture  regime.  Compositions  between  Zr/H  ratios  90/10  and 
65/33  reveals  a  ferroelectric  fenoelectiic  transitioa  between  rhombohedral  space  groups  R3m 
R3c.  This  transition  involves  the  oxygen  octahedra  tilt  about  the  <1 1 1>  directions.^  The  aim  of 
this  investigation  is  to  study  the  inter-relationship  between  octahedral  tilt  domain  structures  and  the 
high  temperature  ferroelectric  domain  structiires.  There  have  been  virtually  no  studies  regarding 
the  domain  structures  of  octahedral  tilt  systems  in  perovskites.^  The  transmission  electron 
microscope  offers  an  attractive  meapf  to  study  this  subtle  phentmiena  owing  to  electron  scattering 
factors  being  much  larger  than  x-ray  scattering  factors.  The  diffraction  contrast  imaging  also 
allows  a  direct  means  ID  study  the  domain  states.^ 
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Ceramic  samples  of  Pb(Zr,Ti)03  were  made  using  conventional  solid-state  sintering 
techniques.  The  starting  raw  material  powders  were  of  analytical  grade  quality  and  included  PbO, 
Zr02,  and  Ti02.  Two  compositions  were  made  for  this  study  Pb(Zro.9Tio.i)03  and 
Pb(Zro.6sTio.3S)03-  These  compositions  were  batched  according  to  stoichiometric  ratios  and 
accounting  for  loss  of  ignitions.  The  raw  powders  were  ball  milled  with  ethanol  solvent  for  48 
hours  for  complete  mixing.  Perovskite  phases  were  fully  formed  |||Br  calcining  for  four  hours  at 
900*C,  as  determined  by  a  Sdnag  X-ray  diffrKtometer.  The  calcined  ^wder  was  ball  milled  for 
24  hours.  Binder  and  1  wt%  excess  PbC  was  added  to  80  meshiR^ed  powders.  Green  pellets 
with  60%  theoretical  density  were  formed  using  uniaxial  pressure  followed  by  binder  burnout  at 


550*C  for  1  hour.  Sintering  was  undertaken  at  t250*C  for  2  hours  to  give  pellets  from  91-94% 
th  '  aid  dl 


TEM  observations  were  made  420  SfEM  U 


tilt  stage  was  used  to  make  in  dtu  TEM  observatioos  beiweiil 


|o  a  thickness  -  SO  pm. 
binute  epoxy.  Ion  beam 
I  a  beam  current  -  1  mA. 
iatan  liquid  nitrogen  two- 


Results 

Figures  2  (a)  and  (b)  show  the  [110]  zone  axis  diffiactioa  patterns  revealing  {h  +  l/2,k  + 
1/2, 1  -»■  1/2]  pseudo  cubic  superiatdee  reflections  in  Pb<Zto.9Tio.i)03  and  Pb(Zrt).6sTio.35)03 
respectively.  It  is  found  that  by  heating  the  Pb^Snp.as'Ho.^ljPs  saiiq)le  to  80*0  die  superlattice 
reflection  disappears,  Figure  2^j^|pd  cons^sely  with  the  Gatan 

liquid  nitrogen  stage  the  intemiQ^  ^acnamii^  These  lefieeSeM  Aijlbfore,  associated  with  a 

displacive  phase  transitioa.  The  (h  *  1/2,  k  *  l/2»  1  *  1/2}  difftyiigat  are  ctmsisient  with  the 
neutron  dififlactioQ  study  petfonned  by  Glazer.^  Glazer  predicied  dic|Hgin  of  this  superstructure 
to  be  oxygen  octahedral  dlts  within  die  sinqile  petovsldte  structure.  The  oxygen  octahedral  shifts 
with  equal  components  about  the  pseudo-q^ic  perovskite  axis  as  to  give  an  effective  clockwise 
and  anticlockwise  rotatkm  of  oxygen  octahedra  about  the  <in>  directions  parallel  to  the 
ferroelectric  dqtole  diqilaoeiiients  of  die  R3m  idiase,  Rgure  X 

Figure  4(a)  showi  die  daik  field  digEractioo  contrast  image  associated  with  the  superiattice 
reflectioo  in  a  Pb(Zio.9Tio.i)03  subgndn.  A  dark  ribbon-like  boundary  is  observed  under  diese 


imaging  conditions,  and  this  is  believed  to  be  a  wall  separating  out-of-phase  octahedral  tilt  variants. 
Figure  4(b)  shows  the  same  subgrain  imaged  under  a  multiple  beam  bright  field  condition.  This 
reveals  180*  domain  boundaries  within  this  region.  There  is  only  a  weak  spatial  perturbation 
between  the  180*  ferroelectric  domain  walls  and  the  octahedral  antiphase  boundaries. 

Figure  S(a)  shows  a  bright  field  image  which  reveals  ferroelectric  twin  structures  and 
inversion  180*  ferroelectric  domains  typical  in  a  rhombohedral  ferroelectrics.^  Figure  S(b)  shows 
the  same  crystallite  imaged  in  dark  field  with  a  superlattice  reflecdon.  The  antiphase  boundary 
contrast  is  again  observed  and  shows  a  strong  interaction  with  the  twin  boundaries.  Generally  we 
found  antiphase  boundaries  are  terminated  on  the  twin  boundaries,  grain  boundaries,  or 
alternatively  contained  within  closed  loops.  The  antiphase  boundaries  in  the  Pb(Zr.Ti)03  ceramics 
predominantly  terminate  on  either  twin  boundaries  or  grain  boundaries.  Region  X  in  Figure  5(b) 
shows  an  example  of  the  antiphase  boundary  tominated  at  a  { 1 10)  twin  domain  boundary  region 
Y  shows  an  antiphase  boundary  to  be  coincident  with  a  (001 )  domain  wall. 

Discussion 

From  the  above  results  we  can  infer  that  the  R3m  ferroelectric  phase  has  only  twin  and 
inversion  domains.  Twin  domains  being  109*  or  71*  type  and  twin  on  habit  planes  { 1 10)  and 
(100),  respectively.  These  observations  are  consistent  with  earlier  observations  on  modified 
rhombohedral  Pb(Zr,Ti)03  ceramics.'^  The  octahedral  tilt  transition  is  driven  by  a  zone  boundary 
soft  mode  resulting  in  the  doubling  of  the  unit  cell.  The  reflections  observed  ate  consistent  with 
the  proposed  model  by  Glaaer,  Figure  3.  This  transition  gives  rise  to  two  additional  variants 
which  are  separate  from  each  other  with  antiphase  boundaries.  The  antiphase  boundaries  ate 
slightly  perturbed  by  180*  domain  walls  and  are  strongly  perturbed  by  the  twin  boundaries.  In 
some  incidents  we  noted  that  die  anti-phase  boundaries  were  consistent  with  the  twin  walls  which 
may  imply  a  coqiliiig  between  the  gradients  of  tlw  tilt  and  the  ptdarization. 
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Figure  Captions 


Figure  1.  The  PbZrOs-PbTiOa  Phase  Diagram  [JafTe,  Cook,  Jaffe  (1971)]. 

Figure  2.  [110]  Zonei^  Axis  Diffraction  Patterns  (a)  Pb(Zro.9Tio.i)03  and 

(b)  Pb(Zro.6sTio.3s)03,  both  at  room  temperature,  and 

(c)  Pb(Zra.i^.3S)03  at  SOX. 

Figure  3.  Schematic  representation  of  the  oxygen  octahedral«tiiting  in  R3c 
phase. 

Figure  4  (a)  Dark  field  image  of  octahedral*tilt  superlattice  reflection  revealing 

anti^phase  boundaries  (APB);  (b)  multiple  bright  field  image  of  same 
area  revealing  typical  180*  domain  wall  contrast. 

Figure  S.  (a)  Dark  field  image  of  ferroelectric  domain  walls;  (b)  Dark  field 
image  of  superlattice  reflection  revealing  interaction  of  APBs  witK 
twin  donudn  walls. 
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TEMPERATURE  (*C) 


Figure  1.  The  PbZrOs'PbTiOs  Phase  Diagram  [iaffe*  Cook»  Jaffa  (1971)]. 
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Figure  2.  1110]  Zone  Axis  Diffraction  Patterns  (a)  Pb(Zrs.fTie.i)03  and 

(b)  Pb(Zrt.esTie.35)03,  both  at  room  temperature,  and 

(c>  Pb(Zr«.(5cTfn  at  liO'C. 


Figure  3.  Schematic  representation  of  the  oxygen  octabedral*tilting  in  R3c 
phase. 


Figure  4  (a)  Dark  field  image  of  octahedraMilt  superlattice  reflection  revealing 

anti'phase  boundaries  (APB);  (b)  multiple  bright  field  image  of  same 
area  revealing  typical  180*  domain  wall  contrast. 


Figure  5.  (a)  Dark  field  image  of  ferroelectric  domain  walls;  (b)  Dark  field 

image  of  superlattice  reflection  revealing  interaction  of  APBs  with 
twin  domain  walls. 


